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ABSTRACT 

We present the properties of the first 250 ftm blind sample of nearby galax¬ 
ies (15 < D < 46Mpc) containing 42 objects from the Herschel Astrophys- 
ical Terahertz Large Area Survey (i7-ATLAS). HerscheVs sensitivity probes 
the faint end of the dust luminosity function for the first time, spanning a 
range of stellar mass (7.4 < < 11.3 logic Mq), star formation activity 

(-11.8 < SSFR < —8.9 logic yr gas fraction (3-96 per cent), and colour 

(0.6 < FUV-Ar^ < 7.0 mag). The median cold dust temperature is 14.6 K, colder than 
in the Herschel Reference Survey (18.5 K) and Planck Early Release Compact Source 
Catalogue (17.7 K). The mean dust-to-stellar mass ratio in our sample is higher 
than these surveys by factors of 3.7 and 1.8, with a dust mass volume density of 
(3.7 ± 0.7) X IO^MqMpc”^. Counter-intuitively, we find that the more dust rich a 
galaxy, the lower its UV attenuation. Over half of our dust-selected sample are very 
blue in FW-Ks colour, with irregular and/or highly fiocculent morphology; these 
galaxies account for only 6 per cent of the sample’s stellar mass but contain over 
35 per cent of the dust mass. They are the most actively star forming galaxies in the 
sample, with the highest gas fractions and lowest UV attenuation. They also appear 
to be in an early stage of converting their gas into stars, providing valuable insights 
into the chemical evolution of young galaxies. 


Key words: galaxies: general - galaxies: irregular - galaxies: evolution - galaxies: 
ISM - submillimetre: galaxies - infrared: galaxies 


1 INTRODUCTION 
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On average, half of all starlight emitted by galaxies is ab¬ 
sorbed by dust and thermally re-emitted in the Far-InfraRed 
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(FIR) and submillimetre (submm) (Fixsen et aL, 1996; 
Driver et al., 2007). Dust is particularly prevalent in star 
forming regions, where the high-energy photons emitted by 
young stars are highly susceptible to absorption by dust 
grains (Fitzpatrick, 2004). The thermal emission from dust 
in galaxies is normally dominated by the hot component, 
which is mostly heated by star-forming regions (Kennicutt, 
1998; Kennicutt et ah, 2009). Thermal emission from dust 
therefore provides an invaluable avenue for the study of star 
formation. The cold diffuse dust component dominates the 
mass of dust in galaxies (Draine et ah, 2007; Law et ah, 2011; 
Ford et ah, 2013; Hughes et ah, 2014), but it is unclear if 
this cold component is also indirectly heated by star forma¬ 
tion through Ultraviolet (UV) photons leaking from birth 
clouds, or if the evolved stellar population is mainly respon¬ 
sible (Bendo et ah, 2012; Boquien et ah, 2011; Bendo et ah, 
2014). Ultimately, the ratio of recent/evolved stellar heating 
is likely to depend on an individual galaxy’s dust geometry 
and star formation activity (Dunne, 2013). Knowledge of 
how this ratio depends on measurable properties (e.g. mor¬ 
phological type, Ltir, colour, etc) would make the determi¬ 
nation of star formation rates from FIR measurements more 
reliable. 

The Interstellar Medium (ISM) is enriched by evolved 
stars, which synthesise heavy elements and then introduce 
them to the galactic environment. Interstellar dust is now 
understood to be the product of both winds from evolved 
stars (Ferrarotti V Gail, 2006; Sargent et ah, 2010), and of 
core-collapse supernovae (SNe), the end-point in the fleeting 
lives of massive stars (Dunne et ah, 2003, 2009; Barlow et ah, 
2010; Matsuura et ah, 2011; Gomez et ah, 2012b; Indebe- 
touw et ah, 2014). However studies of both local (Matsuura 
et ah, 2009; Dunne et ah, 2011) and high-redshift (Morgan 
V Edmunds, 2003; Dwek et ah, 2007; Michalowski et ah, 
2010; Rowlands et ah, 2014b) galaxies have shown a dis¬ 
parity between the rate at which dust is removed from the 
ISM (either by star formation or interstellar destruction), 
and the rate at which stars replenish it. As such, the origin 
of dust in galaxies is still very much an open question. 

It is difficult to develop a thorough understanding of 
galaxies without also understanding the properties of their 
ISM. As FIR and submm astronomy has matured, numer¬ 
ous projects have been undertaken to characterise dust in 
galaxies. The galaxy dust mass function was first mea¬ 
sured for ~ 200 InfraRed (IR) and optically selected galax¬ 
ies by the Submillimetre Gommon-User Bolometer Array 
(SGUBA) Local Universe Galaxy Survey (SLUGS, Dunne 
et ah, 2000; Vlahakis et ah, 2005). This is being followed in 
the era of the Herschel Space Observatory^ (Pilbratt et ah, 
2010) by the Herschel Reference Survey (HRS, Boselli et ah, 
2010) and the Key Insights in Nearby Galaxies Far-Infrared 
Survey with Herschel (KINGFISH, Kennicutt et ah, 2011). 
However, these and other FIR surveys of nearby galaxies 
may have been hindered by the fact they are not dust se¬ 
lected, instead they are selected for their properties at other 
wavelengths. Large-area missions such as with the Infrared 
Astronomical Satellite (IRAS, Neugebauer et al. 1984) and 


^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 


more recently Planck (Planck Gollaboration et ah, 2011a) 
provide blindly selected EIR/submm samples of local galax¬ 
ies, including the recent sample by Glemens et al. (2013), 
but lack resolution and sensitivity when compared to the 
targeted surveys. 

Now, however, with the advent of blind, large-area sur¬ 
veys such as the Herschel Astrophysical Terahertz Large 
Area Survey (i7-ATLAS, Eales et ah, 2010) we finally have 
an unbiased and unrivalled view of the dusty Universe, with 
resolution and sensitivity hitherto only found in targeted 
dust surveys. 

In this paper, we use H-ATLAS to select local dusty 
galaxies, and investigate the properties of sources chosen on 
the basis of their dust mass. In Section 2 we introduce the 
observations and sample selection. In Section 3 we give an 
account of our extended-source photometry. In Section 4, we 
discuss the key properties of our local i7-ATLAS sample. In 
Section 5 we compare the properties of our sample with other 
samples of nearby dusty galaxies. In Section 6 we examine 
the gas and dust evolution of the galaxies in our sample. 
A companion paper on the dust properties of Hl-selected 
galaxies in the local Universe will be presented in De Vis et 
al. (m prep.). We adopt the cosmology of Planck Gollabo¬ 
ration et al. (2013), specifically Ho = 67.30kms“^ Mpc“^, 

= 0.315, and Qa = 0.685. 

2 HERSCHEL DATA AND THE HAPLESS 

SAMPLE 

2.1 Observations 

Observations for H-ATLAS were carried out in parallel mode 
at 100 and 160 /xm with the Photodetector Array Gamera 
and Spectrometer (PAGS, Poglitsch et ah, 2010) and at 250, 
350 and 500 /xm with the Spectral and Photometric Imag¬ 
ing REceiver (SPIRE, Griffin et ah, 2010) instruments on 
board Herschel. Descriptions of the H-ATLAS data reduc¬ 
tion can be found in Ibar et al. (2010) for PAGS, and Pas- 
cale et al. (2011) and Valiante et al. (m prep.) for SPIRE. 
Photometry in the SPIRE bands was performed upon maps 
reduced for extended-source measurements. Our i7-ATLAS 
PAGS maps were reduced using the Scanamorphos (Rous¬ 
sel, 2013) pipeline, with appropriate corrections made for 
the relative areas of the reference pixels on the focal plane. 

This work makes use of the H-ATLAS Phase-1 Version- 

3 internal data release (Valiante et ah, in prep., Bourne et 
ah, in prep.), which comprises 161.6 deg^ coincident with 
the Galaxy And Mass Assembly (GAMA, Driver et ah, 
2009) redshift survey. GAMA provides spectroscopic red- 
shifts, along with supplementary reductions and mosaics of 
ultraviolet (UV) GALEX (Morrissey et ah, 2007; Liske et 
ah, submitted.; Andrae et ah, in prep.), optical SDSS DR6 
(Adelman-McGarthy et ah, 2008), Near-InfraRed (NIR) 
VISTA VIKING (Sutherland, 2012), and Mid-InfraRed 
(MIR) WISE (Wright et ah, 2010; Gluver et ah, 2014) data; 
details of these reprocessed maps can be found in Driver et 
ah, {in prep.). 


© 2015 RAS, MNRAS 000, ??-?? 
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The source extraction algorithm used in iTATLAS 
(MADX, Maddox et ah, in prep.^ Valiante et ah, in prep.) 
isolates >2.5 a peaks in the SPIRE 250/xm maps and then 
measures the fluxes in all three SPIRE bands at the position 
determined by the 250/xm fit. Eor our catalogue we further 
select only those sources which have a > 5 a detection at 
250 /xm. 

Optical counterparts to if-ATLAS sources were found 
by matching ff-ATLAS sources to SDSS DR7 objects 
(Abazajian et ah, 2009) within a 10'' radius using a likeli¬ 
hood ratio technique (Smith et ah, 2011). This method uses 
the optical-submm separation, SPIRE positional errors, and 
r-band magnitudes of potential counterparts, to derive the 
probability that a given optical galaxy is genuinely associ¬ 
ated with the SPIRE source in question (see Smith et ah, 
2011 and Bourne et ah, in prep, for details of the method). 
Sources with a probability of association R > 0.8 are deemed 
to be ‘reliable’ IDs. 


2.2 The Sample 

A sample of 42 galaxies was assembled from the ff-ATLAS 
Phase-1 Version-3 catalogue in the distance range 15 < D < 
46 Mpc. We wished to sample a volume local enough that we 
retained sensitivity to the lowest-mass and coldest sources, 
populations not previously well studied, and our upper dis¬ 
tance limit of 46 Mpc serves this purpose well. We do not 
include galaxies at D < 15 Mpc, where recessional velocity 
is no longer a reliable indicator of distance. These galaxies 
form the Hersehel-ATLAS Phase-1 Limited-Extent Spatial 
Survey, hereafter referred to as HAPLESS. Multiwavelength 
imagery of the full sample can be found in Appendix A, Eig- 
ure Al. 

We require all sources to have reliable SDSS counter¬ 
parts {R > 0.8, Smith et ah, 2011) and to have been as¬ 
sessed as having science quality redshifts (nQ > 3, Driver 
et ah, 2011) by GAMA. We eyeballed the ff-ATLAS maps 
at the location of all optical sources within the redshift range 
and found no other candidates which may have been missed 
by our ID process. The total number of false IDs expected 
in our sample can be estimated by summing (1 — if) (where 
R is the reliability assigned in the likelihood ratio analysis), 
which gives a false ID rate of 0.7percent. 

Distances were calculated using spectroscopic redshifts, 
velocity corrected by GAMA (Baldry et ah, 2012) to ac¬ 
count for bulk deviations from Hubble flow (Tonry et ah, 
2000). Eor Ho = 67.30 km s“^ Mpc“^, the distance limits we 
impose correspond to a (flow corrected) redshift range of 
0.0035 < ^ < 0.01. Reliable redshift-independent distances 
were used for the two sources for which they were available; 
the distance to UGG 06877 has been determined using sur¬ 
face brightness fluctuations (Tonry et ah, 2001), and the dis¬ 
tance to NGG 5584 is known from measurements of Gepheid 
variables (Riess et ah, 2011). 

Gomparing r-band absolute magnitude (Table 3) to dis¬ 
tance, as shown in the upper panel of Eigure 1, shows that 
there appear to be fewer galaxies at greater distances, where 
larger volumes are being sampled. This is likely to be due 
to large scale structure (Eigure 2), since the percentage cos¬ 
mic variance on the number counts in the volume sampled 
by HAPLESS is ~ 100 per cent (Driver V Robotham, 2010). 
The total number of sources listed in the NASA/IPAG Ex- 
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Distance (Mpc) 


Figure 1. Upper: Absolute r-band magnitude against distance, 
for the 42 galaxies of the HAPLESS sample. The different colours 
and shapes denote whether the galaxy lies in the GAMA09 (or¬ 
ange hexagons), GAMA12 (green circles), or GAMA15 (purple 
pentagons) fields sampled as part of the R-ATLAS Phase 1 data 
release. Lower: The distance distribution of HAPLESS sources in 
the different fields. 


• HAPLESS GAMA09 

• HAPLESS GAMA12 

• HAPLESS GAMA15 



Distance (Mpc) 


Figure 2. Polar plot of the volume sampled by HAPLESS in 
the GAMA09, GAMA12, and GAMA15 fields (right-to-left). The 
positions of the HAPLESS galaxies are shown (same symbols as 
Figure 1). Also shown are all the optical sources with redshifts 
that place them in the volume. The inner distance limit of D = 
15 Mpc is demarked by the dotted black line. 
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tragalactic Database (NED^) in the same volume as our 
sample is 141; we therefore detect 30 per cent of this pop¬ 
ulation. Note that the three if-ATLAS fields (GAMA09, 
GAMA12, and GAMA15; see Figure 1) contain 1, 16, and 25 
HAPLESS sources respectively, representing detection rates 
of 7 per cent, 24 per cent, and 42 per cent. 

We identified the portion of our sample which is lim¬ 
ited by intrinsic 250 /im luminosity; this gives us a volume 
limited sample above L 250 = 8.9 x 10^^ WHz“^ (correspond¬ 
ing to a 250 /xm flux of 35mJy at a distance of 46 Mpc). Of 
the 42 HAPLESS galaxies, 35 would still be detected were 
they located at the furthest distance of the volume sam¬ 
pled. Eollowing the assumptions detailed in Section 4.1, this 
is equivalent to a dust mass limit of 7.4 x 10^ M© for a dust 
temperature of 14.6 K (the average dust temperature of the 
sample, see Section 4.1). The 7 sources fainter than this limit 
are HAPLESS 5, 13, 15, 22, 24, 41, and 42. These objects 
are included when describing the properties of our sample in 
Section 4 and comparing to other surveys in Section 5 but 
are plotted as hollow circles. We correct for the accessible 
volume of these sources when considering dust mass volume 
densities in Section 5.4. 

Einally, UGG 06877 (HAPLESS 1) hosts an AGN 
(Osterbrock & Dahari, 1983), with a significant contri¬ 
bution from non-thermal continuum emission in the UV 
(Markaryan et ah, 1979). This contaminates our star for¬ 
mation rate estimate for this galaxy, rendering it unreliable. 
We therefore omit HAPLESS 1 from discussions of star for¬ 
mation. The key characteristics of the HAPLESS sample, 
such as their common names, redshifts, distances and mor¬ 
phologies, can be found in Table 1. We note that 12 of our 
sources are also part of the smaller nearby sample of H- 
ATLAS galaxies presented in Bourne et al. (2013). 

2.3 Curious Blue Galaxies 

We obtained morphology information from the EEIGI cata¬ 
logue of Baillard et al. (2011), which includes 71 per cent of 
the HAPLESS galaxies; we visually classified the remainder 
(all of which were compact dwarf galaxies) using their pre¬ 
scription. The majority of the galaxies in our sample possess 
very late-type, irregular morphology (Hubble stage T > 8) 
though there are two early types (HAPLESS 1 and 22). Eur- 
thermore, a large fraction of the sample exhibit a high degree 
of flocculence (as defined by the EEIGI catalogue). In all, 
24 of our sample are classed as irregular, and 19 as highly 
flocculent; 31 are one or the other, whilst 11 are both (Ta¬ 
ble 1). These irregular and flocculent galaxies are bright in 
the submm and UV, indicating significant dust mass and 
high specific star formation rates (SSERs). They exhibit ex¬ 
tremely blue UV-NIR colours, arising from the fact that, 
along with being UV-bright, they are NIR-faint; examples 
of this can be seen in Eigure 3. We find a UV-NIR colour-cut 
of FUY-Ks < 3.5 mag to be an effective criterion for identi¬ 
fying such galaxies. This approach is supported by the work 
of Gil de Paz et al. (2007), who found FIJY-Ks colour to be 
a powerful diagnostic for discriminating morphological type. 

These curious blue galaxies with FIJY-Ks< 3.5 span 
a wide range of optical sizes, from 1.3 to 33.3 kpc, with a 


^ http://ned.ipac.caltech.edu/ 



Figure 3. Multiwavelength imagery of four examples of the curi¬ 
ous very blue galaxies found in the HAPLESS sample. From left- 
to-right they are, UGC 09299, NGC 5584, NGC 5733, and NGG 
5705. The bands displayed, from top-to-bottom, are: GALEX 
FUV, SDSS three-colour, VIKING A^-band, and PSF-filtered 
Herschel 250 /xm. Each image is 150" x 150". Note the blue op¬ 
tical colours, flocculent morphologies, NIR faintness, and bright 
extended UV emission. The whole sample is presented in Ap¬ 
pendix A, Figure Al. 


median major axis of 9.3 kpc (derived from r-band R25, the 
radius to the 25^^ magnitude square arcsecond isophote). 
Whilst many of them, particularly the larger examples, pos¬ 
sess disks, they often lack defined spiral structure, and show 
only a weak bulge contribution. 

Whilst the FUV-As colour of UGG 06780 (HAPLESS 
1) is 3.07 mag (which would classify it as a member of the cu¬ 
rious blue population), continuum emission from its AGN is 
contributing to the EUV flux. That said, UGG 06780 clearly 
emits plentiful UV emission not associated with the AGN 
(especially for an early-type), as it is more extended in the 
UV than it is in the optical (Table 3). We therefore opt to 
leave it classed amongst the curious blue population, with 
this caveat. 

GALEX coverage is not available for 2 of the HAPLESS 
galaxies (HAPLESS 19 and 21); however the colour u-Ks is 
well correlated with EUV-A^ (Spearman rank correlation 
coefficient of 0.94 for HAPLESS). By comparing the dis¬ 
tributions of these colours, we can state with 3 a confidence 
that a source with u-Ks< 1.36 will have EUV-As < 3.5. This 
indicates that HAPLESS 19 is a member of our curious blue 
population; visual inspection confirms that it exhibits irreg¬ 
ular and extremely flocculent morphology. 

The EUV-As colours of the HAPLESS galaxies can 
be found in Table 3. Of the 42 HAPLESS galaxies, 27 
(64 per cent) satisfy the very blue EUV-As < 3.5 criterion; 
25 (93 per cent) of these exhibit irregular and/or highly floc¬ 
culent morphology. Of the 15 HAPLESS galaxies with EUV- 
As >3.5, irregular and/or highly flocculent morphology is 
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exhibited by only 7 (47 per cent); a two-sided Fisher test 
suggests this difference is significant at the p < 0.01 level. 

3 EXTENDED-SOURCE PHOTOMETRY AND 
UNCERTAINTIES 

3.1 Extended-Source Photometry 

We conducted our own aperture-matched photometry of the 
HAPLESS galaxies, across the entire UV-to-submm wave¬ 
length range, with exceptions for the IRAS 60 fiva measure¬ 
ments, and for the PACS 100 and 160 /xm aperture fitting; 
these differences are detailed in Sections 3.1.1 and 3.1.2 
respectively. At all other wavelengths, we applied a con¬ 
sistent photometric process, tailored to reliably cope with 
the wide range of sizes and morphologies exhibited by the 
sample across the 20 photometric bands employed. These 
bands are: GALEX EUV and NUV; SDSS ugri, VIKING 
ZYJHKs, WISE 3.4, 4.6, 12 and 22/xm; Herschel-PACS 100, 
and 160/xm; and Herschel-SPIRE 250, 350, and 500/xm. In 
summary, an elliptical aperture was fitted to a given source 
in the EUV-22 /xm bands^. The sizes of these apertures were 
compared to identify the largest, which was subsequently 
then used to perform matched photometry across all bands 
(see Eigure 4). 

In detail, we first cut-out a 2000" x 2000" region centred 
on the target source in each band. In the UV-NIR, bright 
foreground stars were removed. The SDSS DR9 (Ahn et ah, 
2012 ) catalogue was used to identify the locations of the 
brightest ~ 20 per cent of stars in the field. Locations for 
stars in non-SDSS bands were also taken from the SDSS 
catalogue, as it was found to provide the most complete 
and robust identification of the stars present. Each star was 
profiled using a curve-of-growth technique, to determine the 
size of the area to be masked. The pixels in the masked 
region were then replaced by a random sampling of the pixels 
immediately adjacent to the edge of the mask. 

To provide the position angle and axial ratio of the 
source aperture, we identified all of the pixels in the cutout 
that had a SNR > 3 associated with the source, and deter¬ 
mined the vertices of their corresponding convex hull^. As 
the vertices of the convex hull trace the outline of the target, 
least-squares fitting of an ellipse to these points provides the 
position angle and axial ratio (i.e. the shape, but not the size) 
of the elliptical source aperture for the band in question. 

The semi-major axis of the source aperture was de¬ 
termined by placing successive concentric elliptical annuli 
(with the already-determined postion angle and axial ratio) 
on the target, centred on the optical SDSS position, with 
semi-major axes separated by one pixel-width, until a mean 
per-pixel SNR < 2 was reached. As flux associated with a 
source with a Sersic profile will fall beyond the edge of any 
practical SNR cutoff^, the fitted aperture was multiplied by 
a factor of 1.2, large enough to be confident of encompassing 


^ SPIRE bands were not used to define the aperture size due to 
the high levels of confusion noise. 

The convex hull is the tightest polygon that can enclose a given 
set of points. 

^ This is true not only for our SNR technique, but also a curve- 
of-growth approach (Overcast, 2010) and the SDSS Petrosian 
method (Blanton et al., 2001). 



Figure 4. Illustration of the stages of our aperture-fitting pro¬ 
cess, using GALEX EUV imagery of galaxy NGC 5584 (HAP¬ 
LESS 14) as an example. Panel 1 shows the inner 500" x 500" 
portion of the cutout centred upon the target source. Panel 2 
shows all of the pixels in the cutout with SNR > 3. Panel 3 shows 
the significant pixels of the target source, contained within their 
convex hull (red points). Panel 4 shows an ellipse fitted to the 
convex hull; this ellipse provides the position angle and axial ra¬ 
tio of the source aperture. Panel 5 depicts the incremental annuli 
used to establish the semi-major axis at which annular flux falls 
to SNR< 2 (thin concentric lines); 1.2 times this distance is then 
used as the semi-major axis of the source aperture (thick line). 
Panel 6 demonstrates the final source aperture (thick line) and 
sky annulus (thin lines). The apertures at all bands for a given 
sources are then compared to select the largest, which is then 
employed for all bands. 


nearly all the flux, whilst small enough to minimise aperture 
noise. The effects of using different extension factors, tests 
upon simulated sources, and visual inspection, all indicate 
that the factor of 1.2 used here achieves this well. This then 
defined the size of the source aperture. The semi-major and - 
minor axes of the generated apertures were compared across 
wave-bands (after subtracting in quadrature the PSE appro¬ 
priate to that band), and the largest selected as the defini¬ 
tive photometric aperture, to be employed in every band for 
a given source. GALEX EUV or NUV served as the defin¬ 
ing band for most sources, except in the case of early-type 
galaxies (and the more early-type spirals), for which it was 
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generally VIKING Z-band. We also determined the r-band 
R2b and FUV i?28 (the radius to the 25^^ and 28^^ mag¬ 
nitude per square arcsecond isophotes, respectively) of each 
galaxy, by interpolating between the mean surface density 
within annuli of one pixel-width; these values are given in 
Table 3. 

For the FUV-MIR, we subtracted the background using 
a sky annulus with inner and outer semi-major axes of 1.25 
and 1.5 times that of the source aperture. For the PACS and 
SPIRE data we used a larger inner and outer annulus of 1.5 
and 2 times the source aperture, thus ensuring enough pixels 
were sampled to make a valid estimation of the value of the 
background. In both cases, the average background value 
was calculated by taking the iteratively 3cr-clipped mean of 
all pixels within the sky annulus. 

The photometry from the FUV to Ks-hand was cor¬ 
rected for Galactic extinction in line with the GAMA 
method described in Adelman-McGarthy et al. (2008). 

In the case of NGG 5738 (HAPLESS 22), a dwarf lentic¬ 
ular, emission in the submm and UV is confined to a point 
source at the centre of the galaxy, as is often seen in early- 
types (Smith et ah, 2012b). The standard aperture, defined 
by NGG 5738’s much larger optical disc, yields poor-quality 
photometry in the submm bands due to the aperture con¬ 
taining too much background. We therefore opt to utilise 
Herschel point-source photometry in the case of this one ob¬ 
ject. NGG 5738 is unique amongst our sample - in all other 
cases, sources compact in the UV and submm are compact 
across the spectrum. 

3.1.1 IRAS SCANPI Photometry 

Eor IRAS 60 /xm we used the Scan Processing and Integra¬ 
tion Tool (SGANPI®), following the procedure laid out by 
Sanders et al. (2003). The SGANPI tool is unable to process 
non-detections where the estimated background is greater 
than the measured flux; in those cases we record a flux of 0, 
with an uncertainty equal to the IRAS 60 /xm 1 a sensitivity 
limit of 58mJy (Riaz et ah, 2006). 

3.1.2 Herschel PACS Photometry 

In the standard H-ATLAS PAGS 100 and 160 /xm data re¬ 
duction (Valiante et ah, in prep.), Nebuliser (an algorithm 
to remove the background emission, Irwin, 2010) was used to 
flatten the maps after they were run through Scanamorphos 
(which deals with 1// noise on the maps, Roussel, 2013). Eor 
sources with apertures > 2.5^ we used the raw Scanamorphos 
maps instead, as Nebuliser removes some emission at these 
scales. Nonetheless, we still find that using the same aper¬ 
tures for PAGS as for the other bands results in poor pho¬ 
tometry. Elux at 100 and 160 /xm tends to be concentrated 
towards the centres of galaxies, often resulting in a small 
patch of flux at the centre of a much larger aperture; this 
can drive up the aperture noise enough that a source with 
clearly-visible flux can count as a ‘non-detection’. As a re¬ 
sult, we define our PAGS apertures separately, using the 
250 /xm maps for each source, as these are reliable indica¬ 
tors of where dust emission is present. Apart from using a 

® Provided by the NASA/IPAC Infrared Science Archive: http: 
//irsa.ipac.caltech.edu/applications/Scanpi/ 


different band to define the apertures, PAGS photometry 
otherwise proceeds in the same manner as described in the 
main part of Section 3.1. 

3.1.3 Comparison with CALEX-CAMA Photometry 

Given the importance of the UV photometry to this work, 
and the fact that our apertures in most cases were defined 
by analysis of surface photometry in the EUV, we have 
made a detailed comparison of our EUV photometry with 
the Gurve-of-Growth (GoG) EUV photometry provided by 
the GALEX-GAMA survey (Liske et ah, submitted.; Andrae 
et ah, in prep.), which has been extensively used in stud¬ 
ies of GAMA galaxies. The comparison was conducted for 
a subset of 17 HAPLESS galaxies relatively unaffected by 
shredding in the SDSS-based GAMA input catalogue used 
by the automated GALEX-GAMA GoG analysis. Our EUV 
apertures were very similar to those derived by the GALEX- 
GAMA GoG, while our EUV integrated fluxes were initially 
found to be systematically higher by ~10 per cent, with a 
similar degree of scatter. This moderate systematic differ¬ 
ence in integrated flux was traced to differences in approach 
to masking foreground stars in the two methods. The only 
other detectable difference was the additional random un¬ 
certainty (~10 per cent root-mean-square) being introduced 
by our use of Swarped images in place of the individual tiles 
used by GALEX-GAMA. We can conclude that both these 
independent methods are in acceptable agreement. 

3.2 Uncertainties 

To estimate aperture noise for a source, we first 3 cr-clipped 
the pixel values in a given 2000'' x 2000" cutout (exclud¬ 
ing those pixels within the source aperture). Then random 
apertures were placed across the cutout (again excluding the 
location of the source aperture itself). Each random aper¬ 
ture was circular, with the same area as the source aperture, 
and was background-subtracted in the appropriate manner 
for each band, as detailed above. The pixel values in each 
random aperture were inspected; if more than 20 per cent 
lay beyond the cutout’s calculated 3 a threshold, then that 
random aperture was rejected. This process was repeated 
until 100 random apertures had been accepted. We found 
this clipping technique to be necessary in order to prevent 
the final aperture noise estimates being too dependant upon 
the locations of the random apertures; otherwise the pres¬ 
ence of bright background sources in the random apertures 
could cause the aperture noise estimate to vary wildly be¬ 
tween repeat calculations on a given cutout. The WISE 3.4 
and 4.6 /xm maps were found to be particularly vulnerable 
to this effect, due in part to anomalies in the maps (halos, 
etc) caused by bright foreground stars. 

Once 100 random apertures had been accepted, the 
flux in each was recorded, and the standard deviation of all 
100 fluxes was taken to represent the aperture noise. This 
method of aperture noise estimation includes the contribu¬ 
tion from confusion noise in Herschel bands. 

We wanted the uncertainty values of our flux measure¬ 
ments to include not only the background noise and ran¬ 
dom photometric uncertainty, but also include the uncer¬ 
tainty in our ability to measure the total flux of a galaxy. 
To that end, we performed two tests. Eirstly, we repeated 
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the photometry with an aperture size 20 per cent larger for 
each source. Ideally, the fluxes obtained using these larger 
apertures would be identical to those obtained from the nor¬ 
mal apertures; the amount of deviation between the two 
lets us gauge the effectiveness of both our aperture-fitting 
and our background-subtraction. Secondly, we repeated the 
photometry, but instead estimated the background using 
a sigma-clipped median within the sky annulus, instead of 
a sigma-clipped mean. These should both be equally valid 
methods, and so the deviation between the final fluxes re¬ 
turned by them allows us to gauge the limits of our ability 
to accurately determine the background. The additional un¬ 
certainty added by these tests is smaller than the instrumen¬ 
tal calibration uncertainties (see below), except in the opti¬ 
cal bands, where the instrumental calibration uncertainty is 
very small. 

No systematic difference in measured flux was found for 
either of these tests. For each of the two tests, the associated 
error value was determined by calculating the root-median- 
squared deviation across all 42 sources. For each band, these 
two error values were then added in quadrature to the band’s 
calibration uncertainty - as given by Morrissey et al. (2007) 
for GALEX, the SDSS DR9 Data Release Supplement^ for 
SDSS, Edge & Sutherland (2013) for VIKING, the WISE 
All-Sky Data Release Explanatory Supplement for WISE®, 
the PAGS Observers’ Manual^ for PAGS, and the SPIRE 
Observer’s Manual^® for SPIRE (see also Bendo et al. 2013). 
This was then added in quadrature to the aperture noise to 
provide the hnal photometric uncertainty. 

Eor the IRAS 60 /xm photometry acquired separately 
using SGANPI, the reported flux uncertainty is added in 
quadrature to a 20 per cent calibration uncertainty (Sauvage, 
2011 ) to provide the total photometric uncertainty for each 
source. 

The final fluxes and uncertainties in all bands can be 
found in Table Al in Appendix A. 

4 PROPERTIES OF THE HAPLESS GALAXIES 

4.1 Modified Blackbody SED Fitting 

To estimate the dust masses and temperatures of the HAP¬ 
LESS galaxies, we fit Modified BlackBodies (MBBs) of the 
form Su oc v^B(v^Td) to the EIR and submm Spectral En¬ 
ergy Distributions (SEDs), where 13 is the dust emissivity in¬ 
dex. We first tried using a single-temperature MBB, keeping 
/3 fixed at a value of 2 and fitting only those data points with 
A > 100 /xm. This is because the mid-IR part of the SED has 
contributions from very small grains which are transiently 
heated by single photons, and therefore not in equilibrium 
with the radiation field (Boulanger V Perault, 1988; Desert 
et ah, 1990). This contribution results in a power-law be¬ 
haviour for the portion of the SED between 12-70/xm and 
including this data in the single-temperature MBB fit would 


^ http://www.sdss3.org/dr9/ 

® http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/ 

^ http://herschel.esac.esa.int/Docs/PACS/htnil/pacs_om. 
html 

http://herschel.esac.esa.int/Docs/SPIRE/html/spire_ 
om.html 



Figure 5. Example dust SED of UGC 09470 (HAPLESS 30), 
with one- and two-temperature component modified blackbody 
fits attempted (upper and lower panels respectively). Both used 
a fixed (3 = 2. Grey points represent upper limits in the fitting 
routine. This is an example of a galaxy for which a one-component 
dust model systematically underestimates the flux at both 100 
and 500 /xm, whilst overestimating it at 160/xm. 


bias the temperature high. Eigure 5 (upper) shows an exam¬ 
ple of a single-temperature MBB; overall we found that this 
method systematically underestimated the fluxes at 100 and 
500 /xm, whilst overestimating them at 160 /xm. We demon¬ 
strate this using the stacked residuals between the model 
and the data in Eigure 6. 

The residuals suggest that a ‘flatter’ SED, produced ei¬ 
ther by a lower value of (3 or by having dust at a range of 
temperatures (Dunne V Eales, 2001; Shetty et ah, 2009), 
would be more suitable. We next tried leaving as a free 
parameter and found a wide range of /3 values (0-4) could 
adequately fit the HAPLESS sources. Whilst this greatly re¬ 
duced the systematic bias, it did not eliminate it. Kelly et al. 
(2012) recently demonstrated that SED fitting routines 
with a given ‘true’ value of /3, can return a wide range of 
fitted values for (3 (see also Smith et ah, 2013); furthermore 
Galametz et al. (2012a) demonstrated that a variable (3 will 
produce less accurate results than using a fixed value. We 
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Wavelength (/xm) 

Figure 6. The mean residual across the sample between the 
model and the data as a fraction of the uncertainty, x, in each 
band, for the one- and two-temperature modified blackbody fits 
(example in Figure 5). The single MBB approach systematically 
overestimates the fiux at 160/xm whilst underestimating it at 
100 and 500 p,m. The error bars show the uncertainty on the data 
points, defined by cr = N~ 2 . 


therefore use a fixed /3 of 2 in this work, as both observa¬ 
tional (Dunne & Bales, 2001; Clemens et aL, 2013; Smith 
et aL, 2013; Planck Collaboration et aL, 2014) and experi¬ 
mental (Demyk et aL, 2013) evidence suggest values between 
1.8-2.0 are appropriate for nearby galaxies. Using {3 = 2 also 
allows us to easily compare our results to other recent Her¬ 
schel and Planck studies (see Section 5). A single MBB only 
provides a useful approximation if the large grains have a 
narrow range of temperatures (Mattsson et aL, in press)^ 
which appears not to be the case for many galaxies in HAP¬ 
LESS (and other FIR surveys; see Mattsson et aL, in press, 
Bendo et aL, 2014). We therefore opt to use an SED model 
which incorporates two temperature components: 

S. = ^[ M^BG,TA + M,B{u,Tc) ] (1) 

where is the flux at frequency v, is the dust mass ab¬ 
sorption coefficient at frequency u, and Me are the hot 
and cold masses, B(iy,Th) and B(iy,Tc) are each the Planck 
function at frequency u and characteristic dust tempera¬ 
tures Th and Tc, D is the distance to the source. At submm 
wavelengths, the dust absorption coefficient Ki, varies with 
frequency as Kjy (x . 

We performed the two-temperature MBB fitting from 
60-500 /xm; the 22 jam point is used as an upper limit to 
prevent unconstrained warm components from being fit¬ 
ted. A x^"iiiiiiiiiiisation routine was used which incorpo¬ 
rates colour-corrections for filter response function and beam 
area^^. Both temperature components were kept within the 
5-200 K range, but were otherwise entirely free. Note that 

The median colour corrections are 0.957, 0.995, 0.990, 1.000, 
1.004, 0.992 at 60, 100, 160, 250, 350, 500 p,m across our entire 
sample. 


for a galaxy with an SED that is well-fit by a single¬ 
component model, this method is free to assign negligible 
mass to one of the dust components, or fit two identical- 
temperature components. In keeping with other i7-ATLAS 
works, we use a value for the dust absorption coefficient of 
^850 = 0.077m^ kg“^ from James et al. (2002), which we 
extrapolate to other wavelengths using a /3 = 2. 

Using the two-temperature SED fitting, we no longer 
encounter any systematic biases in our model fits to the 
data, as can be seen in the lower panel of Figure 6. Figure 5 
shows an example of both one- and two-temperature fits to 
the SED of HAPLESS 30; the two-temperature fits of all our 
sources are displayed in Figure A2. 

Dust masses^^ and temperatures for the HAPLESS 
galaxies are fisted in Table 2. The temperatures of the cold 
dust components range from 9.2 to 25.6 K, with a me¬ 
dian temperature of 14.6 K. The total dust masses range 
from 2.2 x 10^ to 9.5 x 10^ M©, with a median mass of 
5.6 X 10® M©. Uncertainties in the derived dust masses and 
temperatures were estimated by means of a bootstrapping 
analysis, whereby the fluxes were randomly re-sampled ac¬ 
cording to a Gaussian distribution defined by the flux uncer¬ 
tainties, and a best fit was made to the re-sampled SED; this 
was repeated 1,000 times, and the standard deviation in the 
returned fit parameters was taken to represent their uncer¬ 
tainty. All quoted dust masses are the sum of the cold and 
warm components, though the cold component significantly 
dominates the dust mass budget in most of our galaxies (Ta¬ 
ble 2). 

Some galaxies do have SEDs that would be adequately 
fit by a one-component MBB; in such cases, there is a risk 
that using the two-component model could give rise to a 
spurious low-luminosity cold dust component that would 
yield an artificially large dust mass, and low cold dust tem¬ 
perature. We gauged the potential impact of this effect by 
weighting the dust temperatures, according to: 

^ _ McTc + 

I^ergkted- 

However, this only causes a significant change in tem¬ 
perature for the two galaxies with the lowest values of 
Mc/Myj (HAPLESS 25 and 40). The median T^eighted is only 
0.8 K greater than the median Tc, with no significant differ¬ 
ence to any of the trends with temperature reported in this 
work. It is also important to consider that recent work by 
Bendo et al. (2014) has shown that low-luminosity cold dust 
components are present in some galaxies; in such cases, a 
one-component MBB may be an adequate fit to the data, 
but not reflect the actual nature of the dust in a galaxy. 

It is unclear what relationship the systematic 500 /xm 
excess in our single-temperature MBB fits (Figure 6) bears 
to the submm excess seen by many other authors (Galliano 
et al., 2003; Galametz et al., 2012b; Remy-Ruyer et al., 
2013; Giesla et al., 2014; Gross! et al., 2015) - as we also 
see an excess at 100 /xm, and a deficiency at 160 /xm. The 
two-temperature MBB approach is able to account for all 


The median dust mass in our sample is higher than that in 
the overlapping sample of Bourne et al. (2013), this is due to 
differences in the distances used and the photometry method. 
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Table 2. Dust properties of the HAPLESS galaxies. Dust masses (M) and temperatures (T) were derived using a x^-minimising fit 
to a two-component modified blackbody SED model, given by Equation 1. Uncertainties were determined by means of a bootstrapping 
analysis. 


HAPLESS 

Tc 

(K) 

ATc 

(K) 

Tw 

(K) 

AXto 

(K) 

Mc/Mto 

(log) 

AMc/Mto 

(dex) 

Md 

(logio M 0 ) 

AMd 

(dex) 

Ltir 

(logic L 0 ) 

1 

25.6 

1.9 

59.7 

11.4 

2.1 

0.8 

5.4 

0.1 

9.0 

2 

17.2 

1.6 

67.2 

19.3 

3.6 

1.8 

6.0 

0.2 

8.5 

3 

13.5 

2.4 

27.7 

2.9 

1.3 

0.5 

7.2 

0.2 

9.5 

4 

16.7 

4.9 

32.7 

14.2 

1.4 

1.3 

5.7 

0.4 

8.4 

5 

12.9 

2.6 

52.5 

6.8 

3.1 

1.8 

6.0 

0.3 

8.1 

6 

21.7 

1.0 

64.2 

16.4 

3.2 

1.4 

7.9 

0.1 

10.9 

7 

12.2 

2.8 

23.2 

2.4 

1.1 

0.6 

7.4 

0.2 

9.5 

8 

19.6 

1.3 

58.0 

13.8 

3.0 

1.0 

6.4 

0.1 

9.1 

9 

15.0 

1.6 

44.6 

12.1 

2.9 

1.0 

6.7 

0.2 

8.8 

10 

17.5 

2.6 

30.1 

14.0 

1.2 

1.0 

7.3 

0.1 

10.0 

11 

11.3 

1.4 

17.7 

15.6 

1.5 

2.1 

6.9 

0.2 

8.4 

12 

15.7 

2.4 

30.7 

9.8 

1.5 

1.0 

6.4 

0.2 

8.8 

13 

13.2 

3.1 

50.6 

6.5 

2.9 

1.6 

5.7 

0.3 

7.9 

14 

14.6 

2.4 

26.5 

2.9 

1.1 

0.6 

7.4 

0.1 

9.8 

15 

15.5 

4.3 

62.6 

9.0 

3.4 

2.3 

5.5 

0.5 

7.9 

16 

16.1 

3.2 

30.0 

11.2 

1.4 

1.1 

6.7 

0.2 

9.2 

17 

19.6 

6.1 

56.8 

7.9 

1.9 

1.0 

5.3 

0.6 

8.8 

18 

12.8 

2.5 

29.9 

3.1 

1.5 

0.5 

6.7 

0.2 

9.0 

19 

12.3 

1.6 

28.8 

6.3 

2.2 

1.1 

7.0 

0.2 

8.8 

20 

21.2 

2.4 

57.9 

14.2 

2.7 

1.0 

7.5 

0.1 

10.5 

21 

17.4 

0.9 

34.0 

18.6 

2.5 

1.4 

8.0 

0.1 

10.3 

22 

11.5 

2.1 

34.1 

6.8 

1.9 

0.7 

6.0 

0.4 

8.1 

23 

20.7 

1.9 

58.1 

14.8 

2.8 

1.1 

7.6 

0.1 

10.5 

24 

16.3 

3.9 

16.4 

4.5 

5.7 

1.5 

5.7 

0.3 

8.1 

25 

11.7 

1.1 

21.5 

4.2 

0.2 

2.2 

7.2 

0.1 

9.6 

26 

13.0 

1.8 

27.2 

10.2 

1.7 

1.1 

7.4 

0.2 

9.4 

27 

14.2 

1.6 

64.2 

5.3 

4.3 

1.8 

6.2 

0.2 

8.2 

28 

21.5 

4.1 

36.6 

14.8 

1.3 

1.1 

6.9 

0.1 

10.1 

29 

24.4 

1.5 

66.6 

16.8 

2.8 

1.5 

7.6 

0.1 

10.9 

30 

12.4 

2.4 

30.8 

3.9 

1.8 

0.8 

6.7 

0.3 

8.8 

31 

15.8 

3.3 

28.5 

15.8 

1.7 

1.5 

7.2 

0.2 

9.5 

32 

14.1 

2.9 

28.7 

2.5 

1.0 

0.5 

6.6 

0.2 

9.2 

33 

20.8 

8.2 

42.8 

13.8 

2.1 

3.4 

5.7 

0.8 

8.7 

34 

9.2 

2.7 

24.3 

4.6 

1.9 

1.0 

7.2 

0.5 

8.6 

35 

14.5 

1.5 

55.9 

12.3 

3.4 

0.8 

6.7 

0.2 

8.8 

36 

11.2 

1.1 

61.2 

18.6 

4.1 

1.3 

6.8 

0.2 

8.3 

37 

15.3 

2.8 

27.3 

15.5 

1.6 

1.4 

7.3 

0.2 

9.5 

38 

13.3 

4.0 

23.4 

16.0 

1.5 

2.7 

6.2 

0.6 

8.2 

39 

12.5 

2.9 

27.0 

2.5 

1.2 

0.6 

7.1 

0.3 

9.4 

40 

10.8 

6.3 

24.5 

14.5 

0.5 

2.3 

6.7 

0.6 

9.2 

41 

14.6 

2.5 

46.1 

10.9 

2.3 

0.8 

5.8 

0.3 

8.4 

42 

11.7 

2.5 

11.7 

14.8 

3.2 

2.4 

6.3 

0.4 

7.4 


of our systematic residuals without the need for extremely 
cold (<C 10 K) dust components. 

Total infrared luminosities, Ltir^ from 8-1000//m were 
estimated using the best-fit SEDs and extrapolating below 
60 /xm using a power law to account for the luminosity pro¬ 
duced by the transiently heated small grain population. This 
was done by forcing the SED shape in the mid-IR to a power 
law, anchored to the WISE 22 jim flux (or the WISE 12 jim 
flux if this was not available), and the flux at the peak of 
the best-fit SED (see Ibar et al. 2013 for more details). This 
new SED was then integrated to produce Ltir] note that 
the luminosity using this method was on average 14 per cent 
higher than simply integrating the best-fit MBBs from 60- 
500 /xm. The values determined using this method are in 
good agreement with those determined by De Vis et al. (m 


prep.) derived from performing energy-balance modelling of 
the full UV-submm SED with MAGPHYS (da Cunha et ah, 
2008). The resulting Ltir values are listed in Table 2. 

4.2 Stellar Masses 

To determine the stellar masses of the HAPLESS galax¬ 
ies, we follow the method of Zibetti et al. (2009), which 
assumes a Chabrier (Chabrier, 2003) Initial Mass Eunc- 
tion (IME) and uses z-band luminosity along with a rela¬ 
tionship between stellar mass-to-light ratio and g-i colour. 
This method combines stellar population synthesis models 
(Bruzual, 2007) including dust attenuation and compares 
with a sample of nearby galaxies. Stellar masses arrived at 
by this method have a typical uncertainty of 0.1-0.15 dex 
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Table 3. Miscellaneous measured and derived properties of the 
HAPLESS galaxies. Stellar mass is calculated using Equation 3. 


NO 

B abs 

(Mag) 

R2br 

(arcsec) 

L2^fuv 

(arcsec) 

FUV-ATs 

(mag) 

(logio M©) 

1 

-18.0 

32 

33 

3.07“ 

8.8 

2 

-17.2 

11 

17 

2.03 

8.1 

3 

-19.6 

67 

80 

2.13 

9.2 

4 

-17.8 

21 

14 

3.16 

8.8 

5 

-17.2 

28 

13 

3.58 

8.5 

6 

-22.2 

131 

124 

4.51 

10.8 

7 

-20.1 

124 

125 

2.66 

9.5 

8 

-19.0 

36 

37 

2.41 

9.0 

9 

-18.4 

39 

81 

1.35 

8.6 

10 

-20.6 

89 

10 

4.39 

10.1 

11 

-18.4 

54 

56 

3.74 

8.9 

12 

-17.8 

21 

26 

3.08 

8.6 

13 

-16.3 

13 

8 

3.14 

8.1 

14 

-20.2 

96 

92 

2.72 

9.5 

15 

-17.4 

19 

10 

3.74 

8.6 

16 

-18.8 

43 

36 

4.26 

9.3 

17 

-17.7 

21 

21 

1.55 

8.1 

18 

-18.5 

25 

28 

2.21 

8.7 

19 

-19.1 

102 

- 

<3.5^ 

9.2 

20 

-21.0 

115 

34 

7.00 

10.8 

21 

-22.3 

210 

- 

>3.5^ 

11.3 

22 

-18.9 

24 

34 

7.12 

9.7 

23 

-20.7 

86 

87 

4.96 

10.2 

24 

-16.5 

10 

15 

1.82 

7.6 

25 

-21.2 

97 

51 

5.85 

10.6 

26 

-19.9 

75 

82 

2.39 

9.5 

27 

-17.9 

37 

36 

2.90 

8.6 

28 

-20.6 

68 

39 

3.99 

9.8 

29 

-21.7 

93 

53 

4.55 

10.4 

30 

-18.5 

33 

35 

2.24 

8.8 

31 

-20.1 

65 

74 

2.94 

9.6 

32 

-18.2 

18 

15 

3.60 

8.9 

33 

-17.8 

13 

15 

1.58 

8.3 

34 

-18.7 

35 

15 

1.16 

8.6 

35 

-18.9 

36 

46 

2.78 

9.0 

36 

-17.7 

21 

23 

2.32 

8.4 

37 

-20.2 

61 

56 

4.09 

10.0 

38 

-17.5 

14 

13 

2.70 

8.4 

39 

-19.8 

36 

42 

2.34 

9.3 

40 

-18.8 

32 

29 

2.60 

8.9 

41 

-16.5 

26 

35 

0.64 

7.6 

42 

-15.2 

4 

9 

2.47 

7.4 


Note that UGC 06877 (HAPLESS 1) is an AGN (Osterbrock & 
Dahari, 1983), with a contribution from non-thermal continuum 
emission in the UV (Markaryan et ah, 1979). 


Sources UGG 06780 (HAPLESS 19) and NGG 5746 (HAPLESS 
21) do not have GALEX coverage. We use the u-Kg colour to 
infer whether they belong to the curious blue subset.). 


(Cortese et ah, 2012b) modulo uncertainties in the under¬ 
lying population models. Zibetti et al. (2009) caution that 
their approach may not be appropriate where galaxies have 
very young stellar populations (where i would be overesti¬ 
mated) or significant extinction (where i would be underesti¬ 
mated); ie, sources with obvious dust lanes (only seen in 6 of 
the HAPLESS galaxies, Figure Al). As discussed succinctly 
in Taylor et al. (2011), however, variations in extinction (for 
simple dust geometries), the star formation history, metal- 
licity and age only serve to shift galaxies along the [g — i) 
vs M^jLi relationship, such that uncertainties in these pa¬ 
rameters do not produce large errors in the value of stellar 
mass inferred in this way. 

The full formula we employ to calculate stellar mass is: 

M* = (3) 

where M* is stellar mass and Li is z-band luminosity, both 
in Solar units. Stellar masses are listed in Table 3. 

The stellar masses of the HAPLESS galaxies range 
from 2.6 x 10^ to 2.2 x 10^^ M©, with a median mass of 
9.8 X 10® Mq. The Zibetti et al. (2009) method yields stellar 
masses for our sources in excellent agreement with those pro¬ 
duced by the more sophisticated MAGPHYS tool which has 
the ability to model more extincted or highly star-forming 
systems (De Vis et al., in prep.), and are also in agreement 
with the masses derived by GAMA (Taylor et al., 2011). We 
continue to use the colour method in this work in order to 
compare with other nearby FIR surveys (Section 5). 

4.3 Atomic Gas Masses 

We searched the literature for the highest-resolution 21 cm 
observations available for each of the HAPLESS galaxies. 
We found 15 of our sample have observations in the liter¬ 
ature; the instrument and reference for each can be found 
in Table 4. For the remaining sources, we inspected the Hi 
Parkes All-Sky Survey (HIPASS, Meyer et al., 2004; Zwaan 
et ah, 2004; Wong et ah, 2006) catalogue to find HIPASS 
sources within the full-width half-maximum (FWHM) of the 
Parkes beam (14.3^) centred on the positions of the HAP¬ 
LESS galaxies. To avoid the risk of contamination due to 
confusion, we only accepted matches for which there were 
no other known galaxies within 14.3^ radius on the sky, nor 
within 500kms“^ in velocity. This ensures that the matches 
we accept are isolated in Hi. From HIPASS we identify 16 ad¬ 
ditional 21 cm detections associated with HAPLESS galax¬ 
ies. 

For the 11 sources with neither HIPASS nor literature 
Hi detections available. Hi data for 7 were provided by the 
ALFALFA (Arecibo Legacy Fast ALFA, Giovanelli et ah, 
2005) survey (Haynes, priv. comm.). In total we therefore 
have Hi measurements for 38 (90 per cent) of the objects in 
our sample. 

To calculate our Hi masses, we used the standard pre¬ 
scription: 


Mm = 2.36 X (4) 

where Mhi is the mass of atomic hydrogen in Solar units, 
Sint is the integrated 21cm line flux density inJy km s“^, 
and D is the source distance inMpc. 


© 2015 RAS, MNRAS 000 , ??-?? 






12 C.J.R. Clark etal. 


Table 4. Hi properties of the HAPLESS galaxies. The origin column indicates whether 21 cm data comes from the HIPASS or ALFALFA 
catalogues, or published literature values. The Hi centroid velocity Vr and linewidth W 50 are not available for all sources with 21 cm 
measurements. Hi masses were calculated using Equation 4, and gas fraction is defined by Equation 6. Upper limits on the gas mass 
and gas fraction were derived using Equation 5. 


HAPLESS 

^int 

(Jy km s“^) 

Vr 

(km s-l) 

W ^50 

(km s-i) 

Telescope 

Origin 

Mhi 

(logio M 0 ) 

fHI 

J 9 

1 

1.30 

1146 

78 

GET 91m 

Gourtois et al. (2011) 

8.08 

0.13 

2 

1.39 

1308 

- 

Arecibo 

Salzer (1992) 

8.36 

0.60 

3 

23.70 

1387 

222 

Parkes 

HIPASS 

9.56 

0.67 

4 

0.86 

1439 

140 

VLA-D 

Taylor et al. (1995) 

8.24 

0.22 

5 

0.44 

- 

- 

Arecibo 

Impey et al. (2001) 

7.83 

0.18 

6 

72.00 

1462 

306 

Parkes 

HIPASS 

10.16 

0.19 

7 

60.90 

1539 

243 

Parkes 

HIPASS 

10.03 

0.74 

8 

5.70 

- 

- 

Arecibo 

Sulentic & Arp (1983) 

9.08 

0.52 

9 

46.90 

1537 

198 

Parkes 

HIPASS 

9.94 

0.96 

10 

35.80 

1528 

287 

WRST 

Popping & Eraun (2011) 

9.82 

0.31 

11 

5.90 

1624 

187 

Parkes 

HIPASS 

9.17 

0.62 

12 

3.97 

1560 

176 

Arecibo 

ALFALFA 

8.90 

0.66 

13 

0.38 

1713 

26 

Arecibo 

ALFALFA 

7.87 

0.37 

14 

27.10 

1638 

198 

Parkes 

HIPASS 

9.76 

0.47 

15 

1.08 

1652 

72 

Arecibo 

ALFALFA" 

8.34 

0.35 

16 

4.26 

1673 

205 

Arecibo 

ALFALFA 

8.97 

0.32 

17 

3.50 

1749 

120 

VLA-D 

Taylor et al. (1995) 

8.96 

0.88 

18 

- 

- 

- 

- 

- 

<8.67 

<0.45 

19 

26.90 

1729 

225 

Parkes 

HIPASS 

9.85 

0.82 

20 

43.50 

1736 

431 

GET 300 ft 

Davis & Seaquist (1983) 

9.98 

0.12 

21 

30.70 

1724 

556 

WRST 

Popping & Eraun (2011) 

9.83 

0.03 

22 

- 

- 

- 

- 

- 

<8.70 

<0.09 

23 

25.60 

1748 

294 

Parkes 

HIPASS 

9.78 

0.28 

24 

2.89 

1859 

100 

VLA-D 

Taylor et al. (1995) 

8.93 

0.95 

25 

5.30 

- 

- 

Parkes 

Eottinelli et al. (1990) 

9.08 

0.03 

26 

27.90 

1760 

184 

Parkes 

HIPASS 

9.80 

0.66 

27 

8.40 

1836 

224 

Parkes 

HIPASS 

9.31 

0.83 

28 

5.50 

1878 

150 

Parkes 

HIPASS 

9.16 

0.17 

29 

44.50 

1897 

317 

GET 300 ft 

Davis & Seaquist (1983) 

10.07 

0.28 

30 

3.80 

- 

- 

NED 

NED^ 

9.00 

0.62 

31 

13.30 

1891 

177 

Parkes 

HIPASS 

9.64 

0.51 

32 

1.81 

1916 

113 

Arecibo 

ALFALFA 

8.71 

0.39 

33 

6.10 

1973 

60 

VLA-D 

Taylor et al. (1995) 

9.35 

0.91 

34 

6.50 

2033 

99 

Parkes 

HIPASS 

9.37 

0.86 

35 

3.22 

- 

- 

Arecibo 

Schneider et al. (1990) 

9.11 

0.52 

36 

2.04 

2143 

98 

Arecibo 

ALFALFA 

8.93 

0.77 

37 

- 

- 

- 

- 

- 

<9.03 

<0.09 

38 

1.41 

2433 

127 

Arecibo 

ALFALFA 

8.79 

0.71 

39 

8.80 

2510 

148 

Parkes 

HIPASS 

9.55 

0.60 

40 

3.40 

1622 

148 

Parkes 

HIPASS 

8.87 

0.43 

41 

6.40 

1098 

124 

Parkes 

HIPASS 

8.65 

0.90 

42 

- 

- 

- 

- 

- 

<8.76 

<0.96 


“ Classified by ALFALFA as a low SNR source (SNR = 5). 

^ A 21cm Sint value for HAPLESS 30 (UGC 09470) is available on NED, but no reference is provided. Despite this, the corresponding 
Hi properties of HAPLESS 30 are typical of the HAPLESS sample, thus we opt to include it. 


The Hi properties for each source are listed in Table 4, 
the atomic gas masses range from 6.8 x 10^ to 1.5 x 10^° M©, 
with a median mass of 2.3 x 10^ M©. 

The remaining sources fall below the HIPASS detection 
limit, which typically spans the range 1.6 x 10 ® < Mhi < 
9.8 X 10®) M© for the distance range of our sample (Haynes 
et ah, 2011). We determine a 3 cr upper limit on the Hi mass 
on our undetected sources using the following prescription 
from Stevens et al. (2004): 

Mhi < 2.36 x 10 ®/)^( 3 ( 7 )^ 18 ^ 14 ^ (5) 


where a is the RMS noise in a single channel (0.013 Jy), 
D is the distance in Mpc, the accounts for the number of 
uncorrelated channels (the velocity resolution of HIPASS is 
18kms“^), and W50 is the linewidth measured at 50 per cent 
peak intensity. We use the average value of W 50 observed 
for the HIPASS-detected HAPLESS galaxies to estimate the 
upper limits on the Hi mass (Table 4). 

To quantify how gas-rich a galaxy is, we calculate the 
atomic gas fraction for galaxies with detected Hi masses 
(with upper limits quoted for non-detections); this is defined 
as: 
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HI _ MhI 

“ Mhi + M, ^ ^ 

where provides a lower limit on the fraction of the bary- 
onic mass in the gas phase (as molecular gas is not consid¬ 
ered in this work). 

If there is sufficient optical depth in the line of sight 
for Hi clouds, the Hi fluxes and masses could be under¬ 
estimated due to self-absorption. Bourne et al. (2013) show 
this correction is on average a factor of 1.08 for the over¬ 
lapping sample of galaxies between their sources and HAP¬ 
LESS. As we lack the necessary information to calculate 
the self-absorption for other nearby galaxy surveys (see Sec¬ 
tion 5.1), we do not consider self-absorption here, but note 
that our gas masses, particularly for edge on galaxies, could 
therefore be underestimated by this effect. 

Finally, we do not include the molecular gas compo¬ 
nent in this work due to the lack of uniform measurements 
for this sample. Since the molecular component only tends 
to dominate the gas budget in more-massive, ear her-type 
spirals (Saintonge et ah, 2011), our lack of molecular gas 
information is unlikely to make a substantial difference to 
the interpretation in this work. Using the scaling relations 
for H 2 /H 1 and stellar mass from Bothwell et al. (2014), the 
molecular-to-atomic gas ratios in our sample are predicted 
to be negligible (< 0.1) for all but 10 of our sources, with 
the remaining galaxies having ratios between 0.1-0.7. The 
predicted H 2 /H 1 ratios for our curious blue galaxies range 
from 0.016-0.14 with a median of 0.06 - suggesting using 
the atomic gas only is an appropriate estimate of the total 
gas component for these sources. Note that adding molec¬ 
ular gas would only serve to increase the gas fractions in 
Table 4. The gas masses and gas fractions for the detected 
galaxies in our sample will be discussed in more detail in 
Section 5.6. 

4.4 Star Formation Rates 

To estimate star formation rate (SFR), we use the Hirashita 
et al. (2003) method of combining UV and IR tracers, specif¬ 
ically following Jarrett et al. (2013) to combine GALEX 
FUV and WISE 22 /xm measurements to give the total SFR 
as: 

SFR = SFRfuv + SFR22 (7) 

where SFRfuv is the FUV-derived unobscured SFR (cal¬ 
culated using Equation 8), and SFR22 is the 22/xm-derived 
obscured SFR (calculated using Equation 9). All SFR values 
are in units of M© yr“^. 

UV emission traces unobscured high-mass stars, indi¬ 
cating star formation on timescales of ~ 100 Myr (Kenni- 
cutt, 1998; Calzetti et ah, 2005). For SFRfuv, we use the 
prescription of Buat et al. (2008, 2011): 

SFRfuv = 10 ^'^^ufuvLfuv (8) 

where ufuvLfuv is the uLjj luminosity in the GALEX FUV 
waveband^^ in units of bolometric Solar luminosity. Buat 

^Fuv — 1.987PHz 


et al. (2012) hnd the uncertainty in this relation to be 0.13 
dex. It was calibrated using 656 local galaxies (described in 
Buat et al. 2007) with stellar masses greater than 10^° M©, 
and extends down to SFRs of O.O7M0yr“^; as such it in¬ 
cludes a range of actively star-forming and quiescent sys¬ 
tems. The stellar masses of our sample extend to lower val¬ 
ues than the Buat et al. (2007) sample; however the Buat 
et al. (2007) sample does cover the full luminosity, SSFR, 
and colour range (specihcally NUV-r against FUV-NUV) 
exhibited by the HAPLESS galaxies. Note that their SFR 
prescription assumes a Kroupa (2001) IMF; we convert it to 
the Ghabrier IMF (which we use to derive stellar masses) 
using a correction factor of 0.94. 

MIR emission comes primarily from hot dust, heated 
by short-wavelength photons emitted from newborn stars, 
and traces star formation on time scales < 100 Myr (Galzetti 
et ah, 2005; Kennicutt V Evans, 2012). The WISE 22/xm 
SFR relation of Jarrett et al. (2013) was calibrated by boot¬ 
strapping to the Spitzer 24 /xm SFR relation of Rieke et al. 
(2009), and is given by: 

SFR22 = (1 - r?)10“®'^^V22L22 (9) 

where 77 is the fraction of MIR emission originating from 
dust heated by the evolved stellar population, and U 22 L 22 
is the uLi, luminosity in the WISE 22 /xm waveband^^ in 
units of bolometric Solar luminosity. Rieke et al. (2009) es¬ 
timate the uncertainty in their Spitzer 24 /xm SFR relation 
to be 0.25 dex, and hnd it to be accurate at gauging the 
star formation giving rise to thermal dust emission in IR- 
selected galaxies. Jarrett et al. (2013) hnd the scatter in 
their WISE 22 /xm bootstrap to this relation to be negligi¬ 
ble (~1 percent), thanks to the close similarity between the 
Spitzer 24 /xm and WISE 22 /xm passbands. 

The value of 77 will vary from galaxy to galaxy depend¬ 
ing on its current star formation activity and dust geometry. 
77 may be calibrated independently if other tracers of dust- 
corrected SFR are available, or calculated theoretically; val¬ 
ues in the literature for star forming samples range from 
0.17 < 77 < 0.55 (Buat et ah, 2011; Hao et ah, 2011; Smith 
et ah, 2012a; Kennicutt V Evans, 2012). 

We hrst set 77 = 0.17 to be consistent with Buat et al. 
2011, and compare our total SFRs using Equation 7 with 
those derived from SED modelling using MAGPHYS (da 
Gunha et al. 2008, De Vis et al. in prep.). These two tech¬ 
niques produce SFRs offset by a median factor of 1.42 (see 
Figure 7). The likely cause is that 77 = 0.17 is not an accurate 
measure of the fraction of 22 /xm luminosity powered by the 
older stellar population for our sample, whereas MAGPHYS 
allows this fraction to be determined by the energy balance 
between the UV and FIR for each source individually. There 
may also be differences in the prescriptions for SFRfuv be¬ 
tween Buat et al. (2011) and the stellar population models of 
MAGPHYS (taken from Bruzual V Gharlot, 2003; Bruzual, 
2007). Finally, the offset could be explained if a bias existed 
towards a larger transiently-heated small grain population in 
our sample compared to the Rieke et al. (2009) calibration 
data (indeed there is some evidence that the 22 /xm emis¬ 
sion is not correlated with SFR in some iJ-ATLAS galaxies, 

1^22 = 13.64 THz 
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Figure 7. The SFRs calculated using Equation 7 compared to 
those derived by De Vis et al. {in prep.) by fitting the full UV- 
submm SEDs of our sample using MAGPHYS da Cunha et al. 
(2008). The offset between the two prescriptions is by a factor of 
1.42, as indicated by the dashed line. 


Bourne et al. 2013). Modulo the offset, the correlation be¬ 
tween the two SFR estimates is tight, with the exception of 
4 outliers. The 3 galaxies below the scatter are HAPLESS 
9, 33, and 34; these sources have extremely blue FIJY-Ks 
colours (< 2.0), and SFRs which are significantly dominated 
by the unobscured, UV component. The outlier well above 
the line (HAPLESS 25) is at the extreme red end (in terms 
of FVY-Ks) of our sample, and has roughly equal contribu¬ 
tions from UV and 22 /am emission to its SFR using Equa¬ 
tion 7. The SFR prescriptions therefore appear to disagree 
in these extreme regions of the parameter space, though we 
leave this for a future study (De Vis et ah, in prep). In order 
to compare our sources with other nearby galaxy studies, 
including the HRS (for which we do not have full multi¬ 
wavelength data) and the Planck sample of Clemens et al. 
(2013) which uses MAGPHYS (see Section 5) we therefore 
reduce our SFRs from Equation 7 by a factor of 1.42 to be 
consistent. Note that this rescaling factor is well within the 
usual variation found between different SFR prescriptions. 

Adding in quadrature the uncertainties in the UV (0.13 
dex) and MIR (0.25 dex) relations in Equation 7 yields an 
uncertainty of 0.28 dex in the derived total SFRs (this does 
not include the uncertainty in the FUV and 22 jam luminosi¬ 
ties of individual sources). 

We also calculate the specific star formation rate 
(SSFR), the SFR per stellar mass (Table 5). The calculated 
SFRs range from 0.01 to 7 . 12 M 0 yr“^, with a median SFR 
of O. 18 M 0 yr“^. Derived SSFRs range from 1.6 x 10“^^ to 
1.4 X 10“^yr“^, with a median SSFR of 1.3 x 10“^°yr“^. 


Table 5. Star formation properties of the HAPLESS galaxies. 
GALEX FUV (unobscured) and WISE 22 p,m (obscured) star for¬ 
mation rates, SFRpuv and SFR22, are calculated according to 
Equations 8 and 9. Where both GALEX and WISE data exists, 
we combine this (Equation 7) to yield the total SFR. 


NO 

SFRpuv 

SFR 22 

SFR 

SSFR 

(logio 

Me yr- 

L 

(logio yr“^) 

1 

_a 

-1.2 

- 

- 

2 

-1.3 

- 

- 

- 

3 

-0.4 

-0.8 

-0.2 

-9.5 

4 

-1.4 

-1.9 

-1.3 

-10.1 

5 

-1.9 

-2.3 

-1.7 

-10.2 

6 

-0.0 

0.4 

0.7 

-10.1 

7 

-0.4 

-1.0 

-0.3 

-9.9 

8 

-0.7 

-1.4 

-0.6 

-9.6 

9 

-0.7 

-1.6 

-0.6 

-9.3 

10 

-0.6 

-0.3 

-0.0 

-10.2 

11 

-1.4 

-1.8 

-1.2 

-10.1 

12 

-1.3 

-1.8 

-1.2 

-9.8 

13 

-2.0 

-2.8 

-2.0 

-10.1 

14 

-0.3 

-0.5 

-0.1 

-9.6 

15 

-1.9 

-2.2 

-1.7 

-10.4 

16 

-1.3 

-1.4 

-1.0 

-10.3 

17 

-1.0 

-1.4 

-0.8 

-8.9 

18 

-0.8 

-1.4 

-0.7 

-9.5 

19 

- 

-1.4 

- 

- 

20 

-1.3 

0.0 

0.2 

-10.6 

21 

- 

-0.1 

- 

- 

22 

-2.3 

-2.5 

-2.0 

-11.8 

23 

-0.7 

0.0 

0.2 

-9.9 

24 

-1.4 

-2.0 

-1.3 

-9.0 

25 

-0.9 

-1.0 

-0.6 

-11.2 

26 

-0.4 

-1.2 

-0.3 

-9.8 

27 

-1.4 

-2.0 

-1.3 

-9.9 

28 

-0.6 

-0.2 

0.0 

-9.8 

29 

-0.3 

0.6 

0.8 

-9.6 

30 

-0.9 

-1.6 

-0.8 

-9.6 

31 

-0.4 

-1.1 

-0.3 

-9.9 

32 

-1.4 

-1.4 

-1.0 

-9.9 

33 

-0.9 

-1.8 

-0.8 

-9.2 

34 

-0.5 

-1.9 

-0.5 

-9.1 

35 

-0.9 

-1.8 

-0.9 

-9.9 

36 

-1.3 

- 

- 

- 

37 

-0.7 

-1.1 

-0.5 

-10.6 

38 

-1.3 

-2.0 

-1.2 

-9.6 

39 

-0.5 

-1.2 

-0.4 

-9.8 

40 

-0.8 

-1.4 

-0.7 

-9.7 

41 

-1.4 

-1.6 

-1.1 

-8.8 

42 

-2.3 

- 

- 

- 


“ Note that HAPLESS 1 has contamination from non-thermal con¬ 
tinuum emission in the UV; therefore we do not quote a value for 

SFRpuv 
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5 PROPERTIES IN COMPARISON TO OTHER 

DUST SURVEYS OF NEARBY GALAXIES 

We now compare HAPLESS to other surveys of dust in lo¬ 
cal galaxies. In this section, we consider our entire sample; 
however those galaxies that are not in the luminosity-limited 
subset of HAPLESS are plotted in hgures as hollow circles. 
Table 6 summarises the median properties of each of the 
samples, and the results of K-S tests between them. 

5.1 The Reference Samples 

5.1.1 The Herschel Reference Survey 

With its stated objective to be the ‘benchmark study of 
dust in the nearby universe’, the 323 galaxies of the Her¬ 
schel Reference Survey (HRS, Boselli et al. 2010) have been 
observed with resolution and sensitivity unrivalled by any 
previous EIR survey. The HRS chose As-band brightness 
as its selection criteria, because it suffers least from extinc¬ 
tion and is known to be a good proxy for stellar mass. The 
velocity range of the HRS (1050 <V< 1750 km s“^), with 
corrections made to account for the velocity dispersion of 
the galaxies of the Virgo Cluster, corresponds to a distance 
range of 15 < D < 25 Mpc (whereas the HAPLESS distance 
range is 15 < D < 46 Mpc). 

The apparent magnitude limit of the late type galaxies 
in HRS is Ks <12, which equates to an absolute magnitude 
limit between Ks < —17.43 and Ks < —18.54, depending 
on the distance of the source between the HRS limits^^. 
Erom this we can ascertain that between 4 and 15 of the 
42 HAPLESS galaxies would have been insufficiently lumi¬ 
nous in Ks to have been included in the HRS^^. These faint 
HAPLESS galaxies are low stellar mass systems that tend 
to have very blue EUV-A^ colours; 13 of the missing 15 sat¬ 
isfy our FVY-Ks < 3.5 criterion. Galaxies seen by A-ATLAS 
that are faint in Ks, but nonetheless dusty, represent an or¬ 
thogonal population to the HRS, and reveal selection biases 
imposed on targeted dust surveys that A-ATLAS, with its 
blind sample, is not susceptible to. Another difference be¬ 
tween the samples is that the HRS contains numerous early 
type galaxies, partly due to the stellar mass selection, and 
partly due to the extensive overlap (46 per cent) of the HRS 
sample with the Virgo cluster. 

To allow for a direct comparison of HAPLESS to the 
HRS, we determined dust masses and temperatures for the 
HRS galaxies ourselves, using our own SED-htting method 
(as detailed in Section 4.1) and their published PACS^^ 
(Cortese et ah, 2014), SPIRE (Ciesla et ah, 2012), and WISE 
(Ciesla et ah, 2014) photometry, along with IRAS 60 /xm 
data we acquired using SCANPI in the same manner as 

For early type galaxies, a brighter flux limit of Ks < 8.7 is 
applied. 

Only 3 HAPLESS galaxies overlap with the distance range of 
HRS and of these, only one would have been bright enough for 
the HRS selection. 

We corrected the HRS fluxes to account for a recently-flxed 
error in the Scanamorphos pipeline used to create the HRS PACS 
maps. The published HRS fluxes at 100 and 160 /xm were mul¬ 
tiplied by 1.01 and 0.93 respectively, the average change (with 
scatter ~2 per cent) in extended-source flux in maps produced 
with corrected versions of Scanamorphos. 


for the HAPLESS galaxies (described in Section 3.1.1). We 
likewise calculated Ltir values for the HRS using the same 
method as for HAPLESS. 

We note that our dust masses for the HRS galax¬ 
ies are on average a factor ~ 2.2 lower than in Ciesla 
et al. (2014), consistent with their assumed lower value for 
AX500 = 0.1 m^ kg“C 

Smith et al. (2012b) also find that the submm emis¬ 
sion of two HRS sources, the giant elliptical galaxies M87 
and M84, contain significant contamination from their AGN. 
Therefore, we do not attempt to fit the SEDs of these 
sources. 

Eor the Hi masses of the HRS galaxies, we used the val¬ 
ues published in Boselli et al. (2014). The published stellar 
masses of the HRS (Cortese et ah, 2012b) were calculated in 
the same way as our own. The UV GALEX and optical SDSS 
photometry of the HRS is described in Cortese et al. (2012a), 
whilst their NIR A^-band photometry (Boselli et ah, 2010) 
was acquired from the 2-Micron All-Sky Survey (2MASS, 
Jarrett et al. 2000). To calculate the star formation rates of 
the HRS galaxies, we employed the same technique as for 
the HAPLESS galaxies (Section 4.4), for which we used the 
published HRS WISE and GALEX photometry. As for the 
HAPLESS galaxies, we obtain morphologies for the HRS 
from EEIGI (Baillard et ah, 2011). 

5.1.2 Planck 

Negrello et al. (2013) used the Planck Early Release Com¬ 
pact Source Catalogue (ERCSC) (Planck Collaboration 
et ah, 2011b) to assemble a sample of nearby galaxies. Their 
fiux-limited sample contains 234 dusty galaxies brighter 
than 1.8 Jy at 550 /xm, at distances < 100 Mpc (with the 
vast majority lying at z < 0.01); the authors estimate the 
sample to be 80 per cent complete. Clemens et al. (2013) 
have used this sample to perform a study of the properties 
of nearby dusty galaxies. We hereafter refer to this as the 
Planck C13N13 sample. 

Whilst the Planck-selected sample benefits from being 
blind and all-sky (excepting the galactic plane zone of avoid¬ 
ance), Planck suffers from lower sensitivity and resolution 
compared to Herschel (3.8' in contrast to 18"). Only 3 of 
the HAPLESS galaxies exceed the 1.8 Jy 550 /xm fiux limit 
necessary to feature in the Planck C13N13 sample (and none 
of those are members of the curious blue subset). 

Clemens et al. (2013) also derived dust masses and tem¬ 
peratures for their sources by fitting two-component MBB 
SEDs with yd = 2, which is consistent with our method. Eor 
the Planck C13N13 sample, the authors adopted a value for 
the dust absorption coefficient of /^sso = 0.0383 m^ kg“^, in 
contrast to the value in this work of = 0.077 m^ kg“^. 
As a result, we have divided their dust masses by a factor 
of 2.01 to permit comparison. 

The Planck C13N13 stellar masses and star formation 
rates were estimated using the MAGPHYS multiwavelength 
SED-fitting package (da Cunha et ah, 2008), which produces 
stellar masses which agree exceptionally well with the Zi- 
betti et al. (2009) method we employ (De Vis et ah, in prep.)] 
both methods also assume the Chabrier IME. Hi data were 
available for 220 (94 per cent) of the Planck C13N13 galax¬ 
ies (Clemens, priv. comm.). Once again, we use EEIGI mor¬ 
phologies (Baillard et ah, 2011). 
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Figure 9. The dust and stellar mass properties of the HAPLESS, HRS, and Planck C13N13 galaxies. Left: The distribution of dust 
masses. Centre: The distributions of stellar masses. Note that Clemens et al. (2013) derive the stellar masses for the Planck C13N13 
sample using MAGPHYS, whilst for the HAPLESS and HRS samples we use Equation 3; however the stellar masses produced by both 
methods are in excellent agreement with each other for the HAPLESS sample (De Vis et ah, in prep.). Right: The distributions of 
(ie, specific dust mass). HAPLESS contains a much higher proportion of very dust-rich galaxies than either of the other two samples. 


Table 6. Median parameters derived for the local-volume surveys compared in this work, including the very blue (FUV-AT^ < 3.5) subset 
of the HAPLESS sample. Results from Kolmogorov-Smirnov (K-S) tests between HAPLESS and the HRS and Planek C13N13 surveys 
are also shown, indicating the likelihood of the null hypothesis that two samples are drawn from the same underlying population. 


Sample 

FFY-Ks 

Tc 

Md 

M* 

Md/M^ 

SSFR 

Mhi 

Md/Mni 

fHI 

Mb 


(mag) 

(K) 

(10® Mo) 

(IQSMo) 

(10-3) 

(10-“yr-l) 

(1O«M0) 

(10-3) 


(10® Mo) 

HAPLESS 

2.8 

14.6 

5.3 

1.0 

4.4 

12.9 

14.4“ 

3.9 

0.52 

2.5 

// Very Blue 

2.4 

14.2 

4.8 

0.6 

6.5 

20.7 

12.1“ 

2.7 

0.66 

2.3 

HRS 

4.6 

18.5 

4.6 

4.9 

1.2 

4.1 

8.5“ 

6.2 

0.18 

5.5 

Planek 

- 

17.7 

41.9 

17.4 

2.5 

6.9 

36.4“ 

11.6 

0.17 

22.4 

K-S (HRS) 

10-® 

lO-'^ 

0.15 

10-® 

10-® 

10-5 

0.03 

(M 

1 

o 

1—1 

10-5 

(M 

1 

o 

1—1 

K-S (Planck) 

- 

10-3 

10-13 

10-“ 

10-3 

0.01 

10-3 

o 

7 

O 

1—1 

10-'^ 

10-10 


“ Gas masses are available for 90 per cent of the HAPLESS sample (93 per cent of the very blue subset), 81 per cent of the HRS, and 
90 per cent of the Planek G13N13 sample. 


Whilst almost identical sets of observed and derived 
properties are shared by HAPLESS and the HRS, a more 
limited set of parameters is available for Planck C13N13; as 
a result, not all of the following analyses can include the 
Planck sample. 

5.2 Colour and Magnitude Properties 

As described in Section 2.3, we find FVY-Ks colour to be 
an effective way of identifying the subset of curious blue 
galaxies in our sample, using a colour cut of FIJY-Ks <3.5. 
We find that 64 per cent (27) of the HAPLESS galaxies sat¬ 
isfy this criterion, compared to only 27percent of the HRS 
galaxies with FIJY-Ks colours available. Given that the 
HRS is Ks-band-selected, it is to be expected that its galax¬ 
ies will tend to exhibit redder FIJY-Ks colour. The distri¬ 
butions of FFY-Ks colours for HAPLESS and the HRS are 
shown in the upper panel of Eigure 8. Whilst the HRS more- 
or-less equally samples a wide range of FlJY-Ks colours, 
with a median of 4.6 (Table 6), the blindly-selected HAP¬ 
LESS galaxies tend to occupy a much narrower range of 


colours, with a median of 2.8. The distributions are signifi¬ 
cantly different. 

As demonstrated by Gil de Paz et al. (2007), FFY-Ks 
colour is a strong indicator of morphology, as is also seen 
in the central panel of Eigure 8. The very blue FFY-Ks 
colours of the HAPLESS galaxies indicate that the dust- 
selected universe is dominated by very late type galaxies. 

The lower panel of Eigure 8 is a colour-magnitude plot 
constructed using FFY-Ks colour and K^-band magnitude. 
Both the blue cloud and red sequence can be seen in the 
distribution of the HRS, at (3, -19.5) and (8.5, -22); however 
our HAPLESS sample is skewed towards bluer colours such 
that the bimodality is not visible in this sample; indeed, 
many of the HAPLESS galaxies are in fact bluer than the 
blue cloud peak seen in the HRS distribution. 

5.3 Dust and Stellar Mass 

Eigure 9 compares the dust mass distributions of HAPLESS, 
HRS, and Planck G13N13. The effect of the 1.8 Jy flux limit 
at 550 /am in the Planck G13N13 sample is immediately 
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FUV-/s:s (mag) 

Figure 8. Upper: The distributions of FUY-Ks colour for the 
HAPLESS (blue) and the HRS (red) samples. The galaxies of the 
blind HAPLESS sample tend to be significantly bluer than those 
of the AT^-band selected HRS. Central: Morphology against FUV- 
Ks colour for HAPLESS and the HRS. Lower: Absolute AT^-band 
magnitude against EUV-AT^ colour for HAPLESS and the HRS. 
Hollow circles indicate galaxies that are below the luminosity- 
limit of the sample. 


apparent; only galaxies with high dust masses (and a few 
less massive but very nearby galaxies) were bright enough 
to be included in their sample, which has a median dust 
mass of 4.2 x 10 ^ M©. The HAPLESS and the HRS have 
different selection effects but ultimately have comparable 
median dust masses (Table 6 ). 

The three samples also exhibit notably different distri¬ 
butions in stellar mass (Figure 9). The high flux limit of 
the Planck C13N13 sample naturally biases it towards more 
massive galaxies. HAPLESS spans the broadest range of 
stellar masses, but on average has lower stellar mass systems. 
The median stellar masses of the three samples span over 
an order of magnitude, and the combination of lower stellar 
masses, but moderate-to-high dust masses, means the HAP¬ 
LESS galaxies have the highest median M^/M* ~ 4.4 x 10“^ 
(ie, specific dust mass) out of the three surveys (Figure 9, 
Table 6 ). The very blue subset have an even higher median 


dust-to-stellar mass ratio of 6.5 x 10“^; despite accounting 
for only 6per cent of the stellar mass in the HAPLESS sam¬ 
ple, the curious blue galaxies account for over 35per cent of 
the dust mass. 

5.4 The Dust Mass Volume Density 

We now measure the dust mass function (DMF) and dust 
mass volume density for HAPLESS. In this analysis, we con¬ 
sider all 42 galaxies in HAPLESS. For the 7 sources that 
are fainter than the luminosity complete limit, we estimate 
their accessible volumes using the 1/Vmax method (Schmidt, 
1968), while for the luminosity complete subset the accessi¬ 
ble volume is simply that between the 15-46 Mpc distance 
limits of the sample (1, 540 Mpc^). 

The upper panel of Figure 10 compares HAPLESS to 
the dust mass functions of Dunne et al. (2011), Vlahakis 
et al. (2005), and Clemens et al. (2013). The HAPLESS data 
points have had the appropriate corrections from Rigby et al. 
( 2011 ) applied to account for the statistical effects of flux 
boosting and incompleteness (Section 2 . 2 ). The iAATLAS 
Science Demonstration Phase (SDP) result for 0 < z < 0.1 
from Dunne et al. (2011) (orange line in Figure 10) is based 
on the first 16 square degree field of H-ATLAS. Their dust 
mass function shown here includes a correction factor of 1.4 
for the known under-density of the GAMA09 field at z < 0.1 
relative to the average from SDSS (Driver et al., 2011). The 
Vlahakis et al. (2005) DMF (green line) used submm/IRAS 
colour relations from the SLUGS survey to estimate 850 /xm 
fluxes, and hence dust masses, for all IRAS galaxies in the 
PSCz catalogue (Saunders et ah, 2000). In order to translate 
their IRAS plus 850 /xm flux estimate to a dust mass they 
needed to assume a temperature model for the SED, and 
their cold fit assumes a cold dust temperature of 20 K, which 
seemed reasonable at the time based on submm studies of 
IRAS galaxies by Dunne V Eales (2001). The Planck DMF 
from Clemens et al. (2013) is based on the 550 /xm luminosity 
function from Negrello et al. (2013) and uses the same flux 
limited sample we have described in Section 5.1.2. Table 7 
lists the parameters for the different Schechter functions; we 
have corrected all DMFs to the same value of Kd and the 
same cosmology used here. We note that uncertainties in 
the distance measurements of the different galaxy samples 
could cause considerable scatter in the shape of the DMF, 
particularly at the high end, as demonstrated by Loveday 
et al. (1992). This would result in an observed DMF that is 
effectively a Schechter function convolved with a Gaussian. 
However, as the distance uncertainties vary both within and 
between the samples we compare here, we only present the 
observed mass functions in this work. 

Above Md ~ 10 ^ M 0 , the HAPLESS data points agree 
with the Planck DMF but are higher than those from Dunne 
et al. (2011) and Vlahakis et al. (2005). Galaxies with 
Md > 10^ Mq account for 87 percent of the total HAPLESS 
dust mass. Below this mass, the HAPLESS data points are 
in closer agreement with the Dunne et al. (2011) DMF and 
directly probe to lower dust masses than any of the previ¬ 
ous works. Vlahakis et al. (2005) and Planck G13N13 find 
a steeper faint-end slope than Dunne et al. (2011) and this 
work, but their direct sampling of the faint end is 1-2 orders 
of magnitude less than achieved here. With poor statistics 
in all surveys at the low-mass end, the varying estimates of 
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Figure 10. Upper: The HAPLESS dust mass function (with 
appropriate incompleteness corrections applied) compared with 
those of Planck C13N13, Vlahakis et al. (2005), and the H-ATLAS 
Science Demonstration Phase functions from Dunne et al. (2011). 
Lower: The 250 piva luminosity function of HAPLESS compared 
with the 2 ; < 0.1 H-ATLAS samples from Dunne et al. (2011), 
Guo et al. (2014) and the 2 : < 0.2 sample from HerMES (Vaccari 
et al., 2010). The error bars on points represent Poisson uncer¬ 
tainty. The functions are plotted as thin dotted lines in the regions 
where they are extrapolated. All have been adjusted to our and 
cosmology. 


the slope agree with each other within their 1 <7 uncertain¬ 
ties and so we do not consider these differences worrying at 
present. 

The most significant variation in the dust mass function 
between HAPLESS and Dunne et al. (2011) is the excess of 
HAPLESS galaxies around Ma ~ 10^ M©. This could be due 
to two possible effects: cosmic variance, or incompleteness in 
the Dunne et al. (2011) DMF. 

The volume probed by HAPLESS (and also by the 
other surveys at the faint end) is very small and subject to 
a large uncertainty due to cosmic variance (~ 166 per cent, 
Section 2.2). This effect can be explored by comparing the 
250 /xm luminosity functions since this removes the compli¬ 
cation of relating the 250/xm emission to the mass of dust. 


Table 7. The best-fit Schecter function parameters of the various 
dust mass functions and luminosity functions compared in this 
work. All have been scaled to the same cosmology and value of 
nd we employ. 


Literature Dust Mass Functions 

Reference 

a 

M* 

0* 



(Mo) 

(Mpc“^ dex“^) 

Glemens et al. (2013) 

-1.34 

5.27 X 10^ 

11.0 X 10-3 

Dunne et al. (2011)“ 

-1.01 

4.22 X 10^ 

7.19 X 10-3 

Vlahakis et al. (2005) 

-1.39 

6.49 X 10^ 

2.97 X 10-3 


Literature Luminosity Functions 


Reference 

a 

L* 

(WHz-i) 

V 

(Mpc“^ dex“^) 

Dunne et al. (2011)“ 

-1.14 

1.53 X 10^4 

6.00 X 10-3 

Guo et al. (2014)^ 

-1.06 

1.12 X 10^4 

3.70 X 10-3 

Vaccari et al. (2010)^ 

-1.14 

2.19 X 1024 

4.22 X 10-3 


® Note that this incorporates the 1.42 correction factor ap¬ 
plied by Dunne et al. (2011) to account for under density in the 
GAMA09 field (Driver et al., 2011). 


^ Guo et al. (2014) use a modified Schechter function to fit 
their LF, with an additional parameter of cr = 0.30 (explained in 
Saunders et ah, 1990). 

Vaccari et al. (2010) do not provide the parameters to their 
250 /xm Schechter fit; these values represent our best fit to their 
quoted data points. 


Any differences in the 250 LF will purely be due to variations 
in the space density of sources in the different samples. We 
compare HAPLESS to the luminosity functions of previous 
authors in Table 7 and Figure 10 and find good agreement 
(within errors) with the 0 < 2 ; < 0.1 H-ATLAS luminos¬ 
ity function from Dunne et al. (2011) (from 16deg^ Science 
Demonstration Phase data, scaled by their density correc¬ 
tion 1.4 factor; this is an updated version of the LF pre¬ 
sented in Dye et al. 2010) and from HerMES (over 14.7 deg^ 
at z < 0.2, from Vaccari et al., 2010). The LF derived from 
H-ATLAS Phase-1 data (161.6 deg^) in Guo et al. (2014) 
is lower compared to HAPLESS. This measure has not cor¬ 
rected for the known underdensity of the GAMA09 field and 
also uses a brighter optical magnitude threshold for inclu¬ 
sion of sources than Dunne et al. (2011). It is not clear how 
much of a difference this will make (detailed LFs for the 
full i7-ATLAS Phase 1 will be presented in future work) but 
overall, this comparison indicates that the HAPLESS vol¬ 
ume represents a region of fairly typical 250 /xm luminosity 
density and certainly is not significantly overdense relative 
to the density corrected Dunne et al. (2011) values. 

The fact that we find a greater dust mass volume den¬ 
sity than Dunne et al. (2011), despite having the same 
250/xm luminosity density detection limit, must therefore 
be ascribed to a difference in the average ratio of L 25 o/Md 
in the two samples, such that HAPLESS includes dustier 
objects for a given 250 /xm luminosity threshold. The reason 
behind this lies in the relationship between our selection pa¬ 
rameter, L 250 , and our parameter of interest, Md- Whether 
or not we detect a given mass of dust is strongly dependant 
upon the temperature of that dust. This is illustrated by 
Figure 11, which compares the relation between 250/xm lu- 
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minosity and dust mass for HAPLESS and the HRS. The 
relationship found by Dunne et al. (2011) is shown as a 
dashed line. There is a scatter of ~ 1 dex in this relation¬ 
ship, due to dust temperature. The HAPLESS galaxies have 
more dust mass for a given 250 /im luminosity than both the 
HRS and the Dunne et al. (2011) relation, because they are 
colder on average than the galaxies in those samples^®. The 
issue is therefore that no surveys have a truly ‘dust mass 
limited sample’ but rather, in the case of H-ATLAS, we 
have a 250 /xm luminosity limited sample. Our luminosity 
limit of L250 = 8.9 x 10^^ WHz“^ for the HAPLESS com¬ 
plete sub-sample translates to an approximate dust mass 
limit of 7.4 x 10^ M©, using the average HAPLESS dust 
temperature of 14.6 K. But if we instead use the warmest 
and coldest temperatures in our sample (10-25 K) this limit 
becomes ‘fuzzy’ and ranges from 5 x 10® to 3 X 10® Mq. We 
are assuming that a 250 /xm luminosity limited sample is 
equivalent to a dust mass limited sample, when in reality it 
is not. If we consider the volume accessible to a source with 
Md = 10^ M© at a dust temperature of 14 K, compared to 
that for a source with the same dust mass but at a temper¬ 
ature of 20 K, we find that the warmer source with the same 
dust mass has an accessible volume 8 times greater than 
the colder one. The small area in the Dunne et al. (2011) 
work, combined with this effect, may have resulted in an 
incompleteness to colder galaxies at the median redshift of 
sources in the Md ~ 10^ M© bins. At the time of the Dunne 
et al. (2011) work, it was not expected to find many galaxies 
with such cold dust temperatures. Euture work in measuring 
the dust mass function for the full Phase-1 H-ATLAS area 
(Dunne et ah, in prep.) will address this issue and aim to 
correct for it. 

To determine the dust mass density in the HAPLESS 
volume, we use the combined dust mass of the individual 
sources, with our luminosity-incomplete sources weighted to 
account for the fraction of our volume in which they can be 
detected, according to: 


Pd 


Vtot 


( 10 ) 


where pd is the dust mass density, Vtot is the total sample 
volume, and Vacc is the accessible volume for a given source 
(in the case of the luminosity-complete sources detectable in 
our entire volume, we treat Vacc — Vtot)- The resulting dust 
mass density in the HAPLESS volume is pd — (3.7 ± 0.7) x 
10^ M© Mpc“^, where the uncertainty includes errors on the 
dust masses of the individual sources and poisson statistics, 
but not cosmic variance (in keeping with the errors from 
other estimates discussed below). 

We integrate the Schechter fits to the dust mass func¬ 
tions of Clemens et al. (2013); Dunne et al. (2011) and Vla- 
hakis et al. (2005) down to the average HAPLESS dust mass 
limit of 7.4 x 10^ M© to calculate their values of pd- We ac¬ 
count for the difference in nd in the Clemens et al. (2013) 
work and also note that the units for (jV in their Table 2 are 


Note that the Dunne et al. (2011) best-fit line passes through 
the 19 K isotherm in their scatter, and does so here also; in¬ 
dicating that for a given luminosity and dust temperature, we 
would find the same dust mass. 
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Figure 11. Dust mass against 250/xm luminosity for the HAP¬ 
LESS and the HRS, colour-coded by cold dust temperature. 
Also shown are the median dust temperature (solid) for the 
HAPLESS sources and the relationship in Dunne et al. (2011) 
(dashed). Filled circles show the HAPLESS luminosity complete 
sub-sample. 


actually Mpc“^ for their own fits and not Mpc“^dex“^ as 
is written in their paper. The values they quote for (jV for 
Dunne et al. (2011) and Vlahakis et al. (2005) are however 
in the correct units (Negrello, priv. eomm.). We also scale all 
values to reflect the cosmology used in this work. The cor¬ 
responding values of the local dust mass volume density are 
Pd — (3.2 ± 0.6) X 10^ M© Mpc“^ for Clemens et al. (2013), 
Pd = (1.3 ± 0.2) X 10^M©Mpc“^ for Dunne et al. (2011), 
and Pd = 1.1 X 10^ M© Mpc“^ for Vlahakis et al. (2005). The 
quoted uncertainty is estimated by retaining the fractional 
uncertainty of the integrated value quoted in the original 
works (where applicable). Driver et al. (2007) also estimate 
the dust mass density from their study of the H-band lumi¬ 
nosity function assuming a constant ratio of Ls/Md- They 
derive a value of (2.1 — 2.7 ± 0.8) x 10^M©Mpc“^, after 
accounting for our choice of cosmology and the difference in 
hid used by Driver et al. (2007)^®. 

The HAPLESS value of pd = (3.7±0.7) x 10^ M© Mpc”^ 
is compatible with those of Clemens et al. (2013) and Driver 
et al. (2007), and is significantly higher than those of Dunne 
et al. (2011) and Vlahakis et al. (2005). The dust mass 
density is dominated by sources at and above the knee in 
the Schechter function, where HAPLESS measures a higher 
space density than the Dunne et al. (2011) and Vlahakis 
et al. (2005) surveys. We believe this could be due to in¬ 
completeness in the other works in accounting for the very 
coldest dusty galaxies we see in HAPLESS. 

Extrapolation of the DMEs to zero mass suggest that 
2-8 percent of dust mass in the local volume is in galax- 


Clemens et al. (2013) did not account for the h scaling in 
Driver et al. (2007); also the value of K,d used by Popescu & Tuffs 
(2002) is lower than that used here, not equal to ours as stated 
in Clemens et al. (2013). Therefore the Driver et al. (2007) value 
is lower, not higher, than their estimate, though agrees to within 
their Icr errors. 
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Figure 12. The distribution of SSFRs derived for the HAPLESS, 
HRS, and Planck C13N13 samples. Whilst the HRS and Planck 
C13N13 samples show a broad range of SSFRs, the HAPLESS 
galaxies generally occupy of much narrow range of values, of rel¬ 
atively high SSFRs. 


ies below the approximate HAPLESS mass limit, and thus 
HAPLESS presents a highly complete census (albeit with 
small statistics at present) of the dust content of the very 
local universe. Future work exploiting the full iL-ATLAS sur¬ 
vey, which covers 600 square degrees of sky (compared to the 
161.6 square degrees surveyed in this work), will be able to 
address these matters with far greater statistical power. 

5.5 Links Between Star Formation, Colour, and 
Dust Temperature 

We compare the relative rates of star formation activ¬ 
ity between the samples using specific star formation rate 
(SFR/M*), with distributions shown in Figure 12. The 
HAPLESS galaxies tend towards higher SSFRs, with a me¬ 
dian of 1.3 X 10“^°yr“^. Only 15 per cent (34 per cent) of 
the HRS (Planck C13N13) galaxies exhibit SSFRs greater 
than this. The difference in SSFR distrbutions is statisti¬ 
cally significant (Table 6), again highlighting that the HRS 
stellar mass selection appears to under-sample the regions 
of the parameter space where the blindly-selected galaxies 
are found. 

Figure 13 (top) shows the cold dust temperature distri¬ 
butions of the three samples. The HRS and Planck C13N13 
distributions are similar, with medians of 18.5 and 17.7 K 
respectively (Table 6). The HAPLESS distribution, how¬ 
ever, is strikingly different, with a broad peak in the 10-17 K 
range and a median temperature of 14.6 K; 71 per cent (30) 
of the HAPLESS galaxies have dust temperatures colder 
than both the HRS and Planck C13N13 medians. 

The relationship between cold dust temperature and 
galaxy morphology for the three samples is shown in the 
2^^ panel of Figure 13. A strong correlation is present; the 
dust in later galaxy types tends to be much colder than 
in low metallicity dwarf galaxies (Remy-Ruyer et ah, 2013) 
and in earlier types (Skibba et ah, 2011; Smith et ah, 2012c). 
The HAPLESS galaxies are heavily skewed towards the late 



Figure 13. Cold dust temperatures relations for the HAPLESS, 
HRS, and Planck C13N13 galaxies. From top-to-hottom - 
The distribution of cold dust temperatures. 2^^: The relation be¬ 
tween morphological type and cold dust temperature; HAPLESS 
is heavily skewed towards cold late-type galaxies. Stellar mass 
versus cold dust temperature with the relation from Bourne et al. 
(2012). 4^^: MdiMi, against cold dust temperature. Curves rep¬ 
resent different observing limits in due to the cold dust 

temperature for a given value of M.^. Hollow circles indicate galax¬ 
ies that are beneath the luminosity-limit of the sample. 
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Figure 14. Scaling relations with FUY-Ks colour for HAPLESS 
and the HRS. From top-to-hottom - F^: Stellar mass against 
FUV-AT^ colour. 2'^^: (ie, specific dust mass) against 

FUV-AT^ colour, showing the strong relationship between colour 
and dust-richness. 3'^^: SSFR against FUV-AT^ colour; the two 
are tightly related, with our FUV-ATs' < 3.5 colour criterion cor¬ 
responding to a SSFR ~ 1.1 X 10“^° Cold dust temper¬ 

ature against FUV-Ar 5 colour; HAPLESS shows a preponderance 
of cold blue galaxies, which only make up a small fraction of the 
HRS. 
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Figure 15. Scaling relations with Lpuy/L tir for the HAPLESS 
and HRS. Upper: Cold dust temperature against LpuyjL tir- 
Centre: against Lpuy/L rir’^ counter-intuitively, the 

more dust-rich a galaxy, the larger the proportion of FUV photons 
that go unabsorbed. Lower: FUV-AT^ colour against Lpuy /Lpip] 
correlation is very tight for bluer galaxies, but far weaker for red¬ 
der galaxies. 


type and very cold end of the distribution. The 3^^ panel 
of Figure 13 compares stellar mass with cold dust temper¬ 
ature. Only a weak correlation is seen (the Spearman rank 
coefficients are 0.23, 0.39 and 0.23 for the HAPLESS, HRS 
and Planck C13N13 samples respectively; only the latter two 
samples are statistically significant). We note that Bourne 
et al. (2012) find a correlation between dust temperature and 
stellar mass for blue cloud galaxies using a stacking analysis 
on the H-ATLAS data; this trend is plotted as the dashed 
line. 

The last panel of Figure 13 shows a strong inverse cor¬ 
relation between cold dust temperature and this is 

particularly tight for galaxies with cold dust temperatures 
below ~ 15K. We will explore physical connections between 
stellar heating sources and dust temperature next but first 
we consider whether this trend could be related to selection 
biases. As it is always easier to detect a warm dust source 
at a given mass than a colder one, the lack of galaxies in 
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the upper right of this plot cannot be a selection bias. If 
galaxies existed in this part of parameter space (high dust 
content and warm) we would see them. We interpret this 
as a ‘heating limit’ - there is simply not enough stellar ra¬ 
diation available to heat the dust present above the upper 
temperature envelope. On the other hand, the lack of cold 
galaxies with low Md/M^ (lower left) may well be due to 
the detection bias against low dust mass objects with cold 
temperatures discussed in Section 5.4. Warmer galaxies can 
still be above the flux limit with smaller quantities of dust, 
and there is indeed more scatter to lower M^/M* values 
at higher dust temperatures. The curves in Figure 13 show 
our observing limits (for inclusion in the luminosity-limited 
subset of our sample) of MdjM^ as a function of temper¬ 
ature for different values of M*. This explains the lack of 
sources in the lower left region of this plot and the apparent 
‘tightening’ of the relationship at cold temperatures. 

Comparing our parameters with FUY-Ks colour in¬ 
stead of dust temperature in Figure 14 shows very tight 
scaling relations of colour against stellar mass, and 

SSFR. The 2^^^ panel of Figure 14 demonstrates that bluer 
galaxies are consistently more dust-rich. We see that across 
the 3.5 orders of magnitude of MdjM^ sampled by HAP¬ 
LESS and the HRS, no galaxies are so dusty that extinction 
takes over and FIJY-Ks colours become redder. This is in 
stark contrast to the Dust-Obscured Galaxies (DOGs) and 
SMGs observed at higher redshifts, where dust-richness gives 
rise to severe extinction, resulting in red UV-NIR colours 
(Dey et ah, 2008; Calanog et ah, 2013; Rowlands et ah, 
2014a). 

Despite blue FIJY-Ks colours indicating plentiful ongo¬ 
ing star formation (Figure 14, 3^^^ panel) there is no correla¬ 
tion between ‘blueness’ and dust temperature (Figure 14, 4^^ 
panel). However, many of the curious blue galaxies are found 
to have very cold dust temperatures. A possible explanation 
for this is that a large fraction of their UV luminosity escapes 
unabsorbed by dust. In Figure 15 we examine Lpuv/L tir in 
relation to dust temperature, and FFY-Ks colour. 

LfuvILtir indicates the number UV photons escaping a 
galaxy (unabsorbed), relative to the amount of energy which 
is absorbed by dust and thermally re-emitted in the IR. In 
the case where most dust emission is powered by absorption 
of UV photons rather than optical photons, this is equiv¬ 
alent to a measure of the UV transparency. Overall, both 
the HAPLESS and HRS samples show that the cold dust 
temperature is anti-correlated with Lfuv / Ltir^ suggesting 
that the higher the factor of UV radiation that is absorbed, 
the higher the temperature of the cold dust. The very bluest 
HAPLESS galaxies have the highest values of Lfuv/Ltir, 
and display a range of cold dust temperatures. The central 
panel of Eigure 15 shows that, counter-intuitively, the more 
dust-rich a galaxy is (as defined by M^/M*), the smaller the 
fraction of the UV luminosity that suffers dust absorption. 
The combination of dust-richness and low attenuation leads 
to the very cold dust temperatures in the bluest galaxies. 
This could be due to some physical difference in the grain 
population, leading to more efficient emission and/or less 
efficient absorption in the UV; or due to a difference in the 
dust-star geometry in the bluest sources. This is beyond the 
scope of this work but be will explored further using the 
James et al. (2002) method and radiative transfer modelling 
respectively (Dunne et ah, in prep., De Vis et ah, in prep.). 


The lower panel of Eigure 15 demonstrates a tight corre¬ 
lation between Lfuv ILtir and FFY-Ks colour - except for 
galaxies on the red sequence (FFY-Ks > 6). In these systems 
there is a range of Lfuv / Ltir at the same colour. Recall¬ 
ing that Lfuv I Ltir is really only an attenuation measure if 
most IR luminosity is powered by UV photons (as opposed 
to optical ones), this wide range of values for Lfuv / Ltir on 
the red sequence may indicate that the dust heating in this 
population is not dominated by UV radiation. Dust in early 
type galaxies is often acquired during interactions (Gomez 
et ah, 2010; Smith et ah, 2012c; Rowlands et ah, 2012) which 
may produce a range of dust geometries and therefore a wide 
range of values for Lfuv / Ltir- 

In summary, the bluest FFY-Ks sources exhibit the 
highest SSERs, the highest specific dust masses, the lowest 
UV attenuation and often display very cold temperatures. 


5.5.1 What is Heating the Cold Dust Component? 

Dust heating in galaxies can occur in a variety of ways 
(see Kennicutt V Evans, 2012; Dunne, 2013). Warm dust is 
thought to be associated with star-forming dense molecular 
clouds, with newly formed stars heating the dust to tem¬ 
peratures >30K (Kennicutt, 1998; Kennicutt et ah, 2009; 
Bendo et ah, 2010). Gold dust is usually associated with the 
diffuse ISM (Rowan-Robinson V Grawford, 1989; Boulanger 
et ah, 1996; Lagache et ah, 1998; Tuffs V Popescu, 2005; Bo- 
quien et ah, 2011; Bendo et ah, 2012). Most dust resides in 
this diffuse environment (Dunne V Eales, 2001; Draine et ah, 
2007) where it is heated by the general InterStellar Radia¬ 
tion Eield (ISRE), and is often known as a ‘cirrus’ compo¬ 
nent (Rowan-Robinson V Grawford, 1989). The ISRE may 
be largely composed of photons in the optical produced by 
the old stellar population; in this case, the cold dust luminos¬ 
ity would be powered by the old stellar population and not 
young newly formed stars (Tuffs V Popescu, 2005; Boquien 
et ah, 2011; Bendo et ah, 2012). However, it is also possible 
that high energy UV photons from low optical depth star 
forming regions could ‘leak out’ and therefore contribute 
to heating the diffuse dust component (Law et ah, 2011; 
Popescu et ah, 2011; Glemens et ah, 2013; Hughes et ah, 
2014). The dust heating will also depend on the distribution 
of dust and stars within a galaxy and the optical properties 
of the dust (see Eoyle et ah, 2013). Bendo et al. (2014) re¬ 
cently used a large number of sources from the Very Nearby 
Galaxy Survey, HRS, and KINGEISH to show that the rela¬ 
tive contributions of young and evolved stars to dust heating 
varies greatly among nearby spiral galaxies. In this section, 
we wish to investigate the relative importance of both the 
young and the old stellar population in heating the bulk 
dust mass (ie, the cold component). 

Our choice of ‘heating’ parameter is influenced by the 
study of Eoyle et al. (2013) who proposed that the amount 
of star formation (or alternatively old stellar luminosity) per 
unit dust mass should determine the temperature of the bulk 
dust component, not simply the amount of star formation 
and/or old stars (ie, SER/stellar mass) or even their surface 
density. If there is more dust to be heated by a particular 
radiation field then its average temperature will be lower. 
To explore this, Eigure 16 (left column) compares SER/M^ 
(a proxy for energy in star formation per unit dust mass) 
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Figure 16. The influence of star formation and the older stellar population upon the temperature of the cold dust in the HAPLESS 
and HRS galaxies. Upper left: SFR/Md. Lower left: L^g/Md. Upper right: SFR surface density, Tisfr- Lower right: Stellar mass surface 
density, /i*. 


and LKg/Md (a proxy for the energy in old stellar photons 
per unit dust mass) with the cold dust temperature. 

In the upper panel we see that higher values of SFRjMd 
correlate with higher values of Tc with Spearman r coeffi¬ 
cient of 0.74 for the combined surveys (r = 0.71 for HAP¬ 
LESS, r = 0.75 for HRS). This was also seen in Planck 
C13N13. In the lower panel, a positive correlation is seen 
between Lks/^s. and Tc (Spearman coefficient r = 0.69 
for the combined surveys - r = 0.64 and r = 0.67 for the 
HAPLESS and HRS respectively) with most sources clus¬ 
tered together with a well-defined maximum described by 
Tmax ~ The galaxies which scatter well above 

this tight cluster of sources are early-types (E and SO). These 
relationships suggest that both star formation and the old 
stellar population are important contributors to the heating 
of the diffuse dust component in these samples, (first noted 
for the HRS by Boselli et al. 2012). On average, we also see 
that for a given value of Tc, the SFR/Md is higher in the 
HAPLESS sources compared to the HRS, whereas this is 
not the case when comparing 

Kirkpatrick et al. (2014) make the same comparison in 
a study of a sample of resolved star forming spirals with 
average M* = 8.2 x 10^ Mq from the KINGFISH survey. 
They use 500 jam luminosity as a proxy for dust mass, 3.6 
/am luminosity to trace the the older stellar population, and 
Ha+ 24 jam emission to trace SFR in a study of star-forming 
spirals. They find a similar relationship between SFR per 
unit dust mass and temperature; however in contrast to this 
work they find no significant correlation between the old 
stellar luminosity per dust mass and temperature. A key 
difference in our approach compared to theirs is that they 
consider resolved regions within their galaxies and so pho¬ 
tons are required to be absorbed within the same pixel they 
were emitted in order to produce a correlation. It is not clear 



Figure 17. The relationship between L^g and SFR, colour- 
coded by FUY-Ks colour (indicating morphology), for the HAP¬ 
LESS and the HRS galaxies. Hollow circles show galaxies beneath 
the luminosity limit of the sample. 


whether the Kirkpatrick et al. (2014) sample would produce 
the same trends as we see if only the global integrated values 
were considered. 

A sufficiently tight, linear correlation between Lkq 
SFR could make it appear that cold dust temperature is cor¬ 
related with the heating parameter SFR/Md, even if Lrs 
alone was driving the cold dust temperature, with no con¬ 
tribution from star formation (and vice-a-versa). To test for 
this, we plot Lks against SFR in Figure 17. A tight, linear 
relation in this plot could give rise to spurious correlations 
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with temperature in Figure 16. However, the scatter in this 
plot is in fact very large - almost 2 orders of magnitude in 
SFR are possible for a given value of Lks ? with distinct se¬ 
quences of ETGs and LTGs (as indicated by their FIJY-Ks 
colour) visible. Given how weak the correlation in this plot 
is (Spearman rank correlation coefficient of 0.28 with both 
samples taken together), it does not seem possible that it 
could be artificially driving the tight relations in the left- 
hand panels of Figure 16 (which have Spearman rank cor¬ 
relation coefficients of 0.74 and 0.69 respectively). This was 
corroborated by using Monte-Garlo simulations of the Lrs 
against SFR relation to generate ‘spurious’ versions of the 
heating relations in the left half of Figure 16; the ‘spurious’ 
simulated plots were never able to replicate the degree to 
which the actual heating relations are tighter than the Lkq 
vs SFR relation. 

We also show the relations between the surface densities 
of star formation (J^sfr) and stellar mass (/x*) against Tc in 
Figure 16 (right column)^°. The surface density of SFR and 
the old stellar population should be a first approximation to 
the average ISRF contributed by both populations (the en¬ 
ergy per unit area), although this assumption is complicated 
in the case of the SFR surface density, as some fraction of the 
UV radiation will be absorbed locally by dust in the birth 
clouds and contribute to heating the warm dust component 
rather than the cold. 

Both the HAPLESS and HRS (see also Boselli et al. 
2012) show a clear correlation between /x* and Tc {r = 0.44 
and r = 0.63 respectively, with r = 0.65 when the two sam¬ 
ples are combined); the stellar mass surface density is also 
higher for HRS galaxies at a given dust temperature. Eurther 
evidence for heating by the old stellar population comes from 
the sample of iTATLAS galaxies in Bourne et al. (2013) (of 
which 12 overlap with our sources) using independent mea¬ 
surements of correlations with dense and diffuse gas compo¬ 
nents. 

Additionally, the HAPLESS galaxies show a weak, but 
significant correlation between star formation surface den¬ 
sity and cold dust temperature (r = 0.32) while, in contrast, 
the HRS shows no significant correlation (r = 0.04). More 
specifically, the HAPLESS sample shows a range of Tsfr at 
the coldest dust temperatures, but requires a higher Tsfr to 
reach higher temperatures; the HRS sample instead shows 
a range of T>sfr at all temperatures. 

This suggests that while both the young and old stel¬ 
lar radiation fields play a role in heating the dust in both 
samples, the HRS dust heating is more strongly influenced 
by the old stellar population while HAPLESS sources are 
on average more strongly heated by the young stellar popu¬ 
lation. As the average stellar mass and SSER of HAPLESS 
are lower and higher respectively than HRS, finding an ISRE 
dominated by young stars is not surprising in the HAPLESS 
systems. 

5.6 Gas Properties 

Here we compare the gas properties of the three samples. 
Eigure 18 shows the Hi mass distribution of the three sam- 

Surface densities were estimated using the r-band R25 (Ta¬ 
ble 3) to determine the optical radius in kpc, assuming that each 
galaxy is circular as a first approximation. 


pies; 90 per cent (38) of the HAPLESS galaxies, 81 per cent 
(263) of the HRS, and 94 per cent (220) of the Planck 
G13N13 galaxies have Hi data available. Interestingly, the 
median HAPLESS Hi mass of 1.4 x 10^ M© is greater than 
the HRS median - despite the median HAPLESS stellar 
mass being 4 times lower than that of the HRS. Once again, 
the bias of the Planck G13N13 sample towards more massive 
objects is manifest (Table 6). 

The Hi gas fractions (Equation 6) of the HAPLESS 
galaxies have a median value of 0.52 (Table 6), and show 
a relatively flat distribution from 0.03 to 0.96, spanning a 
much wider range than those in HRS or Planck G13N13. Of 
the HAPLESS galaxies with Hi detections, 58 per cent (18) 
have baryonic masses which are in fact dominated by their 
atomic gas component. This is without any consideration 
of molecular gas, the inclusion of which would only serve 
to drive up the gas fractions still further. In contrast, the 
HRS and Planck G13N13 distributions are strongly skewed 
towards lower gas fractions, with medians of 0.18 and 0.17 
respectively. A K-S test suggests that HAPLESS galaxies 
are drawn from a different underlying population in terms 
of gas fraction (Table 6). 

The right hand panel of Eigure 18 shows the baryonic 
masses of the three samples, where Mb — Mhi + M*. This 
measure of galaxy mass may be more appropriate for com¬ 
paring samples where stars make up only a small fraction 
of the total baryonic mass of some of the galaxies. Whilst 
HAPLESS and the HRS have very different distributions of 
stellar mass and Hi mass (Eigures 9 and 18), the differences 
are far less pronounced once we consider baryonic mass (see 
Table 6). ThePlanck G13N13 sample is again limited by its 
high 550 /xm flux limit, primarily sampling galaxies with high 
baryonic masses. In the local Universe - where the largest 
halos have already completed more of their star formation 
- this tends to populate the Planck G13N13 sample with a 
relatively high fraction of passive, high stellar mass and low 
gas fraction galaxies. 

In Eigure 19, we consider the properties of the galax¬ 
ies in relation to their atomic gas richness] ie, Mhi/M^,. 
The top panel shows that bluer FVY-Ks colour is strongly 
correlated with higher levels of gas richness; we also note 
that the edges of this distribution appear to be quite sharp; 
for a given FFY-Ks colour, only a small range of Mhi/M i, 
seems permissible. This is also seen in the correlation be¬ 
tween FFY-Ks colour and SSER (Eigure 14), and together 
with the middle panel, shows that SSER is strongly related 
to gas fraction in the local Universe. This relation will be 
explored further in a companion paper (De Vis et ah, in 
prep.) The curious blue HAPLESS galaxies {FFY-Ks < 3.5) 
are the most gas-rich of all; 80 per cent of those with Hi de¬ 
tections contain a greater mass of Hi than of stars, and their 
median gas fraction is 0.66. Gonversely, all but one of the 
HAPLESS sources with Mhi/M i, > 1 are in the curious blue 
category. It transpires that the {FFY-Ks < 3.5) colour cri¬ 
terion we adopted to identify the curious very blue galaxies 
corresponds to the divide between galaxies whose baryonic 
mass is gas-dominated, and those which are star-dominated. 

The dust-to-gas ratio of the samples (Md/Mni) are 
compared in the lower panel of Eigure 19. Until now we 
have described the HAPLESS galaxies, especially the curi¬ 
ous blue subset, as being very dust-rich compared to other 
EIR surveys, in light of their high values of Md/Mi,. But 
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Figure 18. The atomic gas properties of the HAPLESS, HRS, and Planck C13N13 galaxies. Left: The Hi mass distribution. Centre: 
The gas fraction (Equation 6). The HAPLESS sources have higher gas fractions than seen in other FIR surveys of local galaxies. Right: 
The baryonic mass distribution. Whilst the stellar mass and Hi mass distributions of HAPLESS and the HRS are very different, their 
baryonic mass distributions are rather more similar. 


the HAPLESS galaxies are in fact dust-poor relative to their 
gas mass. The median value (see Table 6) of MdjMm for 
the HRS and Planck C13N13 galaxies are 6.2 x 10“^ and 
1.2 X 10“^ (ie, gas-to-dust ratios of ~ 160 and ~ 90), whilst 
the median for the HAPLESS galaxies is3.9xl0“^ (gas-to- 
dust ~ 260). Furthermore, the median Md/Mni of the curi¬ 
ous blue subset is only 2.7xl0~^-a median gas-to-dust ratio 
of ^ 370. In comparing dust-to-gas ratios of high and low 
gas fraction samples in this way, we do need to worry about 
our lack of molecular gas information. The Planck C13N13 
and HRS may have higher H 2 /Hi ratios than the higher gas 
fraction HAPLESS sources (Saintonge et ah, 2011) and so 
the difference in dust-to-gas ratio may be less when this is 
taken into account. The dust properties in relation to gas- 
richness will be explored further in De Vis et al. (m prep). 

Figure 20 compares baryonic mass to Mhi/M^ for the 
HAPLESS, HRS, and Planck C13N13 galaxies. Across all 
three samples, we see a trend where galaxies with large bary¬ 
onic masses tend to have depleted more of their gas than 
smaller objects, though the Planck C13N13 galaxies tend to 
have higher gas-to-stellar mass ratios for a given baryonic 
mass. As the HRS is essentially a stellar-mass-selected sam¬ 
ple, it is biased towards objects that have already converted 
a large fraction of their gas into stars. The high flux limit of 
the Planck C13N13 sample means that it is biased towards 
more massive galaxies; but being selected by dust brightness, 
it nonetheless tends to select the more ISM-rich examples of 
these massive systems. Our blind submm HAPLESS sample 
favours ISM-rich objects and consistently features the most 
gas-rich galaxies of a given baryonic mass. 

6 THE EVOLUTION OF GAS AND DUST 

Here, we will attempt to explain the dust masses and high 
gas fractions of the HAPLESS sources using a chemical and 
dust evolution model to follow the build up of heavy ele¬ 
ments and dust over time as gas is converted into stars. We 
assume a closed box model as the optimistic case for the 


build up of dust (that is, we do not consider inflows and 
outflows of gas) and instead simply follow the gas (and gas 
fraction fg) as it is converted into stars using a star forma¬ 
tion rate ^{t) and an IMF 0(m) (using the Chabrier, 2003 
IMF consistent with Sections 4.2 and 4.4). More details of 
the model can be found in Appendix B; see also Morgan V 
Edmunds (2003) and Rowlands et al. (2014b). 

We assume two possible scenarios for dust formation by 
stars (see Rowlands et ah, 2014b for a more in depth discus¬ 
sion): firstly, where dust is only contributed via the stellar 
winds of evolved Low-to-Intermediate Mass Stars (LIMS); 
and secondly, where dust is contributed via both LIMS and 
SuperNovae (SNe). Whether the majority of dust in galax¬ 
ies is contributed by LIMS or SNe is a long-standing ques¬ 
tion (see the review in Gomez, 2013), though recent results 
(Gall et ah, 2014) suggest not only do SNe form signifi¬ 
cant quantities of dust, but that also these grains are big 
enough (>1 /am) to survive their journey through the harsh 
reverse shock. We use the dust yields from LIMS consistent 
with FIR observations (Ladjal et ah, 2010) and theoretical 
models (Ventura et ah, 2012). For supernova dust yields, we 
use those of Todini V Ferrara (2001), which are consistent 
with the upper range of dust masses observed in histori¬ 
cal SN remnants including the Grab Nebula (Gomez et ah, 
2012b; Owen V Barlow, submitted.), Gassiopeia A (Dunne 
et ah, 2003, 2009; Rho et ah, 2009; Barlow et ah, 2010), and 
SN1987A (Matsuura et ah, 2011; Indebetouw et ah, 2014). 
Type-la SNe are assumed to be negligible contributors to 
the dust budget (Morgan et ah, 2003; Gomez et ah, 2009, 
2012a). Note that we have no dust destruction in our model 
as we want to follow the maximum build up of dust mass at 
a given time^^ (see Rowlands et ah, 2014b). 


We also do not include models for grain growth (Draine, 2009; 
Rowlands et al., 2014b; Mattsson et al., 2014) since this acts to 
counteract the effects of destruction (Dunne et al., 2011; Asano 
et al., 2013). 
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Figure 19. Trends with Mhi (ie, gas richness) for the HAP¬ 
LESS, HRS and Planck C13N13 samples. Upper: FUV-Ks versus 
Bluer colours are strongly associated with higher gas- 
to-stellar mass fractions; The (FUV-AT^ < 3.5) colour-criterion we 
use to define the curious very blue galaxies transpires to corre¬ 
spond to Mhi ~ Af* (vertical dotted line). Centre: Mhi/M^, 
against SSFR. Lower: M^/Mhi against Mhi/M-),. Hollow circles 
indicate galaxies that are beneath the luminosity-limit of the sam¬ 
ple. HIPASS 3 a upper limits (Equation 5) are shown (dotted 
in the case of galaxies not in our 250 /am luminosity-limit sub¬ 
sample.) 


We use four fiducial Star Formation Histories (SFHs), 
shown in Figure B1 and in Table 8, to model the HAPLESS 
galaxies. These SFHs are (i) SFH A - consistent with the 
Milky Way (Yin et ah, 2009); (ii) SFH B - an exponentially 
declining SFR with initial value of = 0.06Moyr“^ 

and a short burst at ~ 1 Gyr; (iii) SFH C - a faster exponen¬ 
tially declining SFR than B, with initial SFR of 2.4 M© yr“^; 
(iv) SFH D - a scaled-up version of SFR C (multiplied by a 
factor of 20) to illustrate the evolution of a galaxy which is 
consuming its gas more rapidly. Using these fiducial SFHs, 
we follow the evolution of the dust mass relative to the bary- 
onic mass as the gas fraction falls. The initial gas mass is 
set to Mg{0) = 4 X 10^° M© for the Milky Way, and for the 
other models we use the observed gas masses and fractions 



Figure 20. Mhi/M i, against the baryonic mass of the HAP¬ 
LESS, HRS and Planek C13N13 samples. The dotted line indi¬ 
cates Mhi/M i, = 1- Symbols as in Figure 19. 



Gas Fraction 

Figure 21. M^/Mb against gas fraction for the three samples. 
Note that the x-axis of this plot goes from a gas fraction of 1 to 
0. The curves show the results from the chemical evolution model 
for different SFHs (Appendix Bl, Table 8) including SFR A - 
consistent with the Milky Way (purple, Yin et ah, 2009); SFR B 
- exponentially declining SFR with initial value of 0.06 M© yr“^ 
and a burst (orange); SFH C - an exponentially declining rate, but 
with higher initial SFR of 2.4M©yr“^ (tourquoise); and finally 
SFR D - a scaled version of SFR C (x 20, tourquoise dashed). The 
dotted purple line is SFH A (MW) with dust from LIMS only. 
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(Table 4) to derive the initial gas masses (these range from 
Mg{0) = 3-5.5 X 10^ Mo). 

The model results are shown in Figure 21. First, we 
compare the dust evolution with dust only from LIMS (dot¬ 
ted line). Second, including dust from LIMS and SNe in 
combination (solid and dashed lines). None of the mod¬ 
els for the former scenario (ie, without SN dust) reach the 
high levels of Md/Ms observed in the HAPLESS, HRS and 
Planck C13N13 samples; this is in line with results from 
other studies, including Morgan & Edmunds (2003), Mat- 
suura et al. (2009), Dunne et al. (2011), Gall et al. (2011), 
and Rowlands et al. (2014b). With dust from both SNe and 
LIMS, SEH models A-C all sit on the same evolutionary 
track in Eigure 21, due to the models with lower star forma¬ 
tion rates (SEHs B & C) than the MW (SEH A) also having 
lower initial gas masses; ie the models lie on the same con¬ 
stant SFR/Mhi tracks. These are in good agreement with 
the HAPLESS galaxies (at high gas fractions) and the HRS 
galaxies at lower gas fractions. The Planck C13N13 galax¬ 
ies (clustered towards lower gas fractions) have somewhat 
elevated Md/Ms compared to the models presented here. 
When we multiply SEH C by a factor of 20 (SEH D) but 
keep the initial gas mass the same as SEH C, the evolution¬ 
ary path is offset, due to the available gas reservoir being 
consumed faster and dust mass reduced due to astration. 

The evolutionary path suggested in Eigure 21 indicates 
that a galaxy’s dust mass will peak when its gas fraction falls 
to ~ 0.5, as predicted in Eales & Edmunds (1996). Therefore 
this is the stage of a galaxy’s development when it is most 
likely to meet the inclusion threshold of a dust-selected sam¬ 
ple such as HAPLESS - the median gas fraction of which 
is indeed 0.5. The stellar-mass selection of the HRS means 
that it is biased towards galaxies where most of the gas has 
already been converted into stars, hence it severely under¬ 
samples the gas-rich portion of this evolutionary path. Sim¬ 
ilarly, the tendency of the Planck C13N13 sample to mainly 
select more massive galaxies means that it too is biased to¬ 
wards systems with low gas fractions. 

7 CONCLUSIONS 

We have presented a dust-selected sample of nearby galaxies 
drawn from the blind H-ATLAS submm survey and intro¬ 
duced a pipeline used to derive photometry, dust masses and 
other properties for the sample. We have also studied corre¬ 
lations between the dust and other properties of the galaxies 
in our sample, and the HRS and Planck samples. We find 
the following results: 

• A typical source seen by Herschel in this blind survey 
has a cold dust temperature of 14.6 K and dust mass of 
5.6 X 10® M©. 

• HAPLESS galaxies have median M^/M* greater by a 
factor of ~3.7 than the galaxies observed as part of the 
Herschel Reference Survey, and a factor of 1.8 than galax¬ 
ies in the Planck Early Release Compact Source Catalogue. 
The median properties of this sample include: (M*) = 9.8 x 
10® M©, (SER) = 0.2M©yr-\ (SSER) = 1.3 x 10“^® yr”^ 
and are amongst the most actively star forming galaxies seen 
in local EIR and submm surveys. 

• This sample contains a high proportion of very blue 
galaxies (defined as FVY-Ks < 3.5). These are generally ir¬ 


regular and/or highly flocculent; such galaxies tend to be 
UV-bright, NIR-faint, dust rich, and low stellar mass, with 
high specific star formation rates. The median dust-to-stellar 
mass ratio of the very blue subset is ~ 3-5 times larger than 
the Planck and HRS samples. Whilst accounting for only 
6 per cent of the stellar mass in our sample, the bluest galax¬ 
ies in our sample contain over 35 per cent of the dust mass. 

• The dust mass volume density of our sample is (3.7 ± 
0.7) X 10® M© Mpc“®, which is higher than some other esti¬ 
mates, but consistent with the value found by Planck. Much 
of this difference seems to arise from the low dust tempera¬ 
tures of the galaxies in our sample, as the 250 /am luminosity 
function of our sample is in good agreement with surveys of 
larger volumes. Note however that our volume suffers from 
a high cosmic variance of ~ 166 per cent. 

• The HAPLESS galaxies are extraordinarily gas rich, 
particularly the very blue sources. Of the 38 HAPLESS 
galaxies detected in Hi, 21 (55 per cent) have atomic gas 
masses greater than or equivalent to their stellar mass. Their 
median gas fraction is 0.52, and 26 per cent have gas frac¬ 
tions > 0.8. The median gas to dust ratios of these sources 
(> 260) is 1.6-3.0 times greater than for the other local sam¬ 
ples of dusty galaxies. 

• The coldest dust seen in the local universe is consis¬ 
tently associated with galaxies that have lots of star forma¬ 
tion relative to their older stellar population. Despite dust 
being so plentiful in these objects, UV photons apparently go 
unabsorbed - giving rise to their very blue FIJY-Ks colours, 
and colder dust temperatures. Comparing the star formation 
and stellar mass surface densities also shows that dust heat¬ 
ing in galaxies selected by HRS is more strongly influenced 
by the old stellar population whereas galaxies selected in 
HAPLESS are more strongly heated by the young stellar 
population. 

• A chemical and dust evolution model confirms these 
galaxies are simply in an earlier stage of converting their 
gas into stars. The bluest galaxies appear to be the most 
immature; they should therefore provide valuable insights 
into the chemical evolution of young galaxies. 

A blind dust-selected sample in the local universe re¬ 
veals very blue, dusty, and gas-rich galaxies. Despite ac¬ 
counting for roughly half of all dusty galaxies, they have 
been severely under-represented in other EIR and submm 
surveys. We suggest the properties of these blue galaxies are 
in line with their ‘immaturity’, and therefore may provide 
useful analogues to very young, high-z galaxies, though we 
note that the interstellar medium in the HAPLESS sources 
is likely to be different. Resolved atomic and molecular gas 
maps of the bluest sources in this sample, combined with 
radiative transfer modelling, 850 /am observations, and inte¬ 
grated optical spectra, should be able to address this and 
test whether these blue-but-dusty galaxies have different 
grain properties or whether the dust is distributed in a 
‘leaky’ geometry (De Vis. et ah, in prep., Dunne et ah, in 
prep., Smith et ah, in prep.). 
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Table 8. The fiducial star formation histories (A-D) and initial gas masses used in this work to model the HAPLESS galaxies. SFH A 
is from Yin et ah, 2009. Also given is the time (tend) where the SFH is truncated to match the present observed SFRs and gas fractions 
of the HAPLESS sample (Section 4). 


SFH 

9(0) 

(Mo) 

V’(O) 

(Moyr-i) 

Burst? 

(Y/N) 

^end 

(Gyr) 

V’(^end) 

(Moyr-1) 

fg (And) 

A (MW) 

4.0 X 10i° 

10 

N 

20 

0.7 

0.11 

B 

3.0 X 10® 

0.06 
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1.35 

0.029 

0.95 

C 

5.5 X 10® 

2.5 
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2.8 

0.5 

0.64 
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5.5 X 10® 

49 

N 

0.1 

48 
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A PROPERTIES OF THE HAPLESS 
GALAXIES 

Multiwavelength imagery of the HAPLESS galaxies can be 
found in Eigure Al. Our UV to EIR photometry of the HAP¬ 
LESS galaxies, with uncertainties, is given in Table Al. Eig¬ 
ure A2 shows the spectral energy distributions of the HAP¬ 
LESS galaxies. 


© 2015 RAS, MNRAS 000 , ??-?? 
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Figure Al. Multiwavelength imagery of each of the HAPLESS galaxies. The bands displayed, from left-to-right, are: GALEX FUV, 
SDSS gri three-colour, VIKING AT^-band, and Herschel 250 /rm. Each cutout is 100" on a side. HAPLESS 19 and 21 do not have GALEX 
coverage. 


© 2015 HAS, MNRAS 000 , ??-?? 
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Figure A1 — continued 
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Wavelength (//m) 


Figure A2. Spectral energy distributions for the HAPLESS sample. The two-temperature modified blackbody fits are shown in red, 
with the contributions from the warm and cold dust components shown by the dashed curves. The grey 22 yam point was treated as an 
upper limit. Sources with dark blue names satisfied the FUY-Ks < 3.5 colour criterion necessary to be counted amongst the curious blue 
sub-population; sources with light blue names did not. 
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Table Al. Photometry of the HAPLESS galaxies, where indicates cases where no observations were available. The semi-major axis is denoted by a, the position angle by and the 
axial ratio (the semi-major axis divided by the semi-minor axis) by ajh. HerschelSVYKEj fluxes were measured using maps reduced for extended-source photometry, but have not been 
colour-corrected. 
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Figure Bl. The star formation histories used to model the HAP¬ 
LESS galaxies. SFH A follows the evolution of the MW (pink, up¬ 
permost solid track, Yin et al. 2009) with initial gas mass Mg(0) = 
4 X 10^® Mq and initial SFR '0(0) = 10 Mq yr“^. SFH B (orange, 
lowermost solid track) has initial gas mass Mg(0) ~ 3 x 10^ Mq 
and initial SFR 0.058 MQyr“^ with a burst at 1 Gyr superim¬ 
posed on top of the exponentially declining rate. SFH C (orange, 
middle solid track) has Mg{0) = 5.5 x 10^ Mq and initial SFR 
2.4 Mq yr“^ and exponentially declines until reaching a gas frac¬ 
tion of fg ~ 0.6 at 2.5 Gyr. SFR D (turquoise, dashed track) is a 
scaled version of SFR C (x20). 


B THE CHEMICAL EVOLUTION MODEL 

Briefly, the equations to follow the evolution of gas and dust 
in the HAPLESS galaxies are: 


Mtot = Mg + M,, (Bl) 

where Mg is the gas mass and M* is the stellar mass. The 
gas mass evolution with time is described by: 

^ = -^Q + e(i). (B2) 

where 0(t) is the rate at which gas is depleted by the SFR, 
and e{t) is the rate at which it is returned as stars die. 

Assuming that mass loss occurs suddenly at the end of 
stellar evolution at time Tm{Tn) (Schaller et ah, 1992), the 
ejected mass, e(t), from stars is: 

nmu 

— [m — mR{m)] — Tm)0(m)(im (B3) 

where ruRt^rn) (from Prantzos et ah, 1993) is the remnant 
mass and is the mass of a star whose age is that of a 
system where a star formed at {t — Tm) has died at Tm- The 
evolution over time of the mass of metals in the ISM, Mz, 
is described by: 


^^ = -Z(i)V-(i) + e.(i) (B4) 

where Z is the fraction of heavy elements by mass in the gas 


phase. The mass of heavy elements ejected by stars at the 
end of their lives is denoted by 

/ mu 

{[m - mR{m)] Z{t - Tm) + mpz) 

X 0(t — Tm)(l>(Tn)dm (B5) 

Yields from stars (mpz) are taken from the theoretical 
models of Maeder (1992) and van den Hoek & Groenewegen 
(1997). The development of the mass of dust with time is 
described by: 

Q?n ?ni^(?n)] Z(t Tm)dold T 'Gfipzdnew^ 

X ^j;{t — Tm)(l>{m)dm — (Md/M^)0(t) (B6) 

where dust is built up from two sources, the fraction of the 
heavy elements that are recycled through star formation and 
ejected in stellar winds (Sold), and the fraction of new ele¬ 
ments freshly synthesised in stars and ejected in both super¬ 
novae and stellar winds (Snew)', the final term describes dust 
removed from the interstellar medium due to astration. 

We use different star formation histories including one 
consistent with the Milky Way, and others consistent with 
galaxies with low star formation rates throughout their evo¬ 
lution (see Figure Bl and Table 8). The initial gas masses 
are derived from the observed properties of the HAPLESS 
galaxies (ie Mg{0) = Mgjfg, Table 4), we truncate the star 
formation histories when they reach the observed gas frac¬ 
tion. E.g. for HAPLESS 3, we start the model with an initial 
gas mass of 5.53 x 10^ Mq, and with the star formation his¬ 
tory SFR C we reach a gas fraction of fg ~ 0.6 at 2.5 Gyrs 
consistent with the observations. Note that the lower SFHs 
are consistent with the current SFRs of the HAPLESS galax¬ 
ies as derived from their UV and MIR fluxes; the SFHs are 
also compatible with the range of SFHs derived from more 
complex multiwavelength modelling of their SEDs (De Vis 
et ah, in prep.). 
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